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Abstract. We present a method to solve initial-boundary value problems for linear and 
integrable nonlinear differential-difference evolution equations. The method is the discrete 
version of the one developed by A. S. Fokas to solve initial-boundary value problems for 
linear and integrable nonlinear partial differential equations via an extension of the inverse 
scattering transform. The method takes advantage of the Lax pair formulation for both linear 
and nonlinear equations, and is based on the simultaneous spectral analysis of both parts 
of the Lax pair. A key role is also played by the global algebraic relation that couples all 
known and unknown boundary values. Even though additional technical complications arise 
in discrete problems compared to continuum ones, we show that a similar approach can also 
solve initial-boundary value problems for linear and integrable nonlinear differential-difference 
equations. We demonstrate the method by solving initial-boundary value problems for the 
discrete analogue of both the linear and the nonlinear Schrodinger equations, comparing the 
solution to those of the con'esponding continuum problems. In the linear case we also explicitly 
discuss Robin-type boundary conditions not solvable by Fourier series. In the nonlinear 
case we also identify the linearizable boundary conditions, we discuss the elimination of the 
unknown boundary datum, we obtain explicitly the hnear and continuum limit of the solution, 
and we write down the soliton solutions. 

7 October 2008 



1. Introduction and outline 

The development of the theory of infinite-dimensional integrable systems was a remarkable 
advance of mathematical physics over the last forty years. One of the key properties of such 
systems is that they can be written as the compatibility condition of an overdetermined linear 
system, called the Lax pair. In turn, the existence of Lax pair is deeply related to many 
other features of these systems. Among them is the inverse scattering transform (1ST), a 
nonlinear analogue of the Fourier transform which can be used to solve the initial value 
problem (IVP). The 1ST was successfully used in the late 1960's and early 1970's to solve 
IVPs on infinite domains or with periodic or quasi-periodic boundary conditions (BCs) for a 
variety of nonlinear partial differential equations (PDEs), differential-difference fully discrete, 
integro-differential equations, etc. (e.g., see Refs. [2, 7, 10, 19] and references therein). 

Following the solution of IVPs, a natural issue was the solution of initial-boundary value 
problems (IBVPs). After some early results [8, 13, 14, 22], however, the issue remained 
essentially open for over twenty years. Recently, renewed interest in the problem has lead 
to a number of developments (e.g., see Refs. [9, 11, 12, 15, 23-30, 39-41] and references 
therein). Particularly important among these is the method developed by A. S. Fokas [23-30]. 
Fokas' method, which is a significant extension of the 1ST, is based on the simultaneous 
spectral analysis of both parts of the Lax pair A crucial role is also played by a relation called 
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global algebraic relation that couples all known and unknown boundary values. Indeed, it is 
the analysis of the global relation that allows one to express the unknown boundary datum 
in terms of known ones plus the initial datum. Importantly, the method also yields a new 
approach to IBVPs for linear PDEs, which allows the solution of new kinds of problems. 

At the same time, the effort to extend the properties of integrable nonlinear PDEs to 
discrete integrable systems has been an ongoing theme in the last thirty years (e.g., see 
Refs. [2, 3, 5, 6, 20, 21, 31, 35, 36, 42] and references therein). The purpose of this work 
is to show that, mutatis mutandis, an approach similar to that for PDEs can also be used 
to solve IBVPs for linear and integrable nonhnear differential-difference equations (DDEs). 
We demonstrate this claim by solving IBVPs for the discrete analogue of the linear and 
nonlinear Schrodinger equations on the natural numbers. Note that the integrable discrete 
nonhnear Schrodinger (IDNLS) equation is an important model since it arises in a number of 
physical and mathematical contexts (e.g., see references in Ref. [6]). 

The outline of this work is the following. In section 2 we solve the IB VP on the natural 
numbers for the discrete linear Schrodinger (DLS) equation, namely the hnear DDE 

qn+i-2qn + qn-\ „ ,11, 
iqn + ^2 =^ (1-1) 

where qn = qn(f) € C, n e N, / s df jdt denotes time derivative and h is the lattice spacing. 
Then, in sections 3 and 4 we consider the IBVP for the integrable nonlinear counterpart 
of (1.1), namely the IDNLS equation or Ablowitz-Ladik (AL) equation [4, 5], 

iq„ + j-^ v\q„\ (q„+i+q„-i) = (1.2) 

(where as usual the cases v = - 1 and v = 1 will be called respectively focusing and defocusing). 
In particular, in section 4 we discuss the eUmination of the unknown boundary datum, the 
linearizable boundary conditions, and we write down the soliton solutions. Finally, in order 
to appreciate the similarities and differences between the method in the discrete versus the 
continuum case, in section 5 we review the solution of IBVPs for the continuum limits of both 
equations, namely the linear and nonlinear Schrodinger equations, and we discuss explicitly 
the correspondence between the method in the discrete case versus the continuum limit. The 
proof of various statements in the text is confined to the Appendix, which also contains a Ust 
of notations and frequently used formulae. 

In both the linear and the nonlinear problem we will require the initial datum to be 
absolutely summable and the boundary datum qo{t) to be smooth, even though the method 
can be formulated under weaker conditions. The constant h can be eliminated from (1.1) 
and (1.2) via the rescalings t' = t/h^ and q'f^it) = hqnit). Thus, for simplicity we will consider 
the rescaled problems throughout (thus effectively setting h=\); however, we will will omit 
the primes except when considering the limit /z — > to recover the solution of the continuum 
cases. The indended meaning should be clear from the context. Also, for brevity we will 
occasionally omit functional dependences when doing so does not cause ambiguity. 



2. Discrete linear Schrodinger equation 

Here we solve the linear problem (1.1), which serves to introduce some of the tools that will 
be used in the nonlinear case. In section 2.1 we derive a Lax pair for (1.1). Then, in section 2.2 
we solve the IVP and in section 2.3 IBVPs via spectral methods. 
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IVP and IBVPfor DLS via Fourier methods. Let us briefly review the solution of the IVP 
and the IB VP via Fourier methods. Doing so we will serve to introduce quantities that will 
also be used later. Consider first the IVP, namely (1.1) with n e Z and with ^„(0) given. 
We require that the initial datum ^„(0) decays rapidly enough as n — > ±00 to belong to ^^(Z), 
the space of sequences {an]nez such that 2nL-oo l^^nl < °°- Introduce the transform pair as 



qik,t)^ Z qn(t)/z"^ 2 e-'^Vn©, (2.1a) 

n=—oo n=—oo 

qnit) = § ^~'q{z,t)dz = ^ ] e'"'q{k,t)dk, {2Ab) 
lzl=l 

where z = e'^, and the contour |z| = 1 is oriented counterclockwise. The transformation A: — > z 
maps A: e R into |z| - 1 and Im ^ into \z\ > 1 (with k - +ioo corresponding respectively to 
z = and z = 00). Use of (2.1) yields the solution of the IVP in Ehrenpreis form as 

q„(t) = (£ z«-ie-'^(^>'§(z,0)dz= ^ re'("'=-^(^>'>§(A:,0)dA:, (2.2) 
\z\=i ^ 

where the hnear dispersion relation is 

aj(z) = 2-(z+l/z) = 2(l-cosA:). (2.3) 

Now consider the IB VP, namely (1.1) with n e N and t E R^, with q„{0) and ^o(f) given. We 
assume q„iO) e ^^(N) and ^^0(0 ^ C(RJ). Introduce the Fourier sine series and its inverse as 

00 1 
q^%, 0=2 qnim/z" -z"), qnit)-^ § (z"-l /z") q^'\z, t) Aziz , 

n=l 4m 1^1^ J 

Use of this pair yields the solution of the IB VP as 

qnit) = §iz"-l /z")/z c-''^^'^' q^'\z, 0) dz - ^ § iz" - I /z")/z c-'''^'^'giz, t) dz , (2.4) 



where 



g{z,t) = {z-llz)]e'^'^'' qoil'W . 





2.1. A Lax pair for the discrete linear Schrodinger equation 

A Lax pair formulation, first discovered for nonlinear PDEs [33], is also possible for linear 
PDEs, and in fact it is the key to solving a wide class of IB VPs [24, 30]. Here we show how a 
Lax pair for the DLS equation (1.1). can be obtained by taking the hnear hmit of the the Lax 
pair of the IDNLS equation (1.2). (As in the continuum limit, an algorithmic way also exists 
to obtain the Lax pair associated to any linear discrete evolution equation. The corresponding 
formalism will be presented elsewhere.) 

It is well-known that the IDNLS (1.2) is a reduction of the Ablowitz-Ladik (AL) 
system (A.4) [5]. A Lax pair for (A.4) is given by the overdetermined linear system (A.5). To 
obtain the linear hmit of (A.5), let Q„ = 0(e), and take 3>n(z,0 = v„(z,0 = ivi,n,V2,n)' to be a 
two-component vector The leading order solution of (A.5) is then v„(z, f) = Z"e''^"'~' ""^'/Zy^^ 
where \o = (fi.o.fi.o)' is an arbitrary constant vector. Choosing v2,o = 1 and keeping terms up 
to 0(e) then yields the following scalar hnear system for viy. 

vi,„^i -zvu = <7„z-«e-*-i/^>''/2 ^ (2.5a) 

Kn-kz- l/zfvi,n = iizqn-qn-l/z)z-''s-'^'-'"^'^\ (2.5b) 
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Enforcing the compatibility of (2.5a) and {2.5b) now yields the discrete hnear Schrodinger 
equation (1.1). To eliminate the dependence on z" from the right-hand side (RHS) of (2.5a), 
we now perform the rescahng z' = and <^„ = ft-i^i(z-i/z) '/^^^ ^_ Dropping primes for 
simplicity, we then obtain the following Lax pair for (1.1): 

(pn+l-Z(pn = qn, ^„ + icoiz)(pn = K^n- qn-l/z) , (2.6) 

where a)(z) is given by (2.3) as before. Indeed, although it may not be obvious at this point, 
the meaning of the variable z in (2.6) coincides exactiy with that of z in (2. 1). 

The rescaling z' = z^ between the linear and the nonlinear problem is the discrete 
analogue of the rescahng k' - 2k in the continuum limit. Such rescaling will reflect on the 
location of the jumps in the Riemann-Hilbert problem (RHP) for the IB VP in the nonlinear 
problem, which will difi"er from the corresponding locations in the Unear problem. 



2.2. IVPfor DLS via spectral analysis of the Lax pair 

We now solve the I VP for (1.1) using spectral methods. Doing so will introduce some of the 
ideas that will be useful for the IB VP and nonlinear case. Making use of the integrating factor 
as in (2.3)], we introduce the modified eigenfunction 

Mz,t)=z-"e^'^'Mz,t), (2.7) 
which satisfies the following modified Lax pair: 

(A„+i-(A„ = e''^(*<7„/z"+i, ^„ = e''^(*/(<?„ -<?„_! /z)/z" (2.8) 

Of course the above linear system is also compatible if q„it) satisfies (1.1). It is then easy to 
define ^^^'^\z, t) as the solutions of (2.6) which vanish as n — > +oo, respectively: 

<t'n\z,t)- "l. qm{t)f-"'-\ (l>f\z,t) = - qMz"-'"-\ (2.9) 

m=—oo m=n 

Note that d)-}\z,f) is analytic as a function of z for |z| < 1 and continuous on |z| = 1, while 
<pn iz,t) is analytic for |z| > 1 and bounded for |z| = 1. The jump conditions obtained by 
evaluating (f)''^'^\z,t) on |zl = 1 then yield a scalar RHP: 4}i^\z,t)-cf>„{z,t) = z''"'§(z,0, where 
qiz,t) is given by (2.1). However, the difference 4>n'' ~<p'n^ solves the homogeneous version 
of (2.6), and hence it depends on n and t only through the factor z" e'"^'-^'. Evaluating (2.9) at 
{n,t) = (0,0) we can then rewrite the jump condition as: 

0!/\z, t) - cf>^^\z, t) = z"-'e-™«'§(z,0) , Izl = 1 . (2. 10) 

Equations (2.9) imply 0^'\O,O = qn-i{t) + 0, and (^i^'(z,0 ^ as z ^ oo. Thus, the RHP 
defined by (2.10) is trivially solved by applying standard Cauchy projectors, namely: 

0„(z,O= § ^d^, (2.11) 

where the contour is oriented counterclockwise, as usual. Then, inserting (2.11) into the LHS 
of the first of (2.6), one obtains the solution of the IVP as (2.2). 

The continuum limit of (2.2) yields the solution of the linear Schrodinger equation. 
Indeed, reinstating the lattice spacing h, the solution of the IVP for the DLS (2.2) is 

1 njh 

qn{t) = — r e'(«*''-"(^^'^§(A:,0)dA:, (2.12a) 

2^ -iih 

where now oj{k) = 2(1 - cos kh)/h^ and 

qik,t)=h i e-'"%(0. (2.12fc) 

n~-oo 

Then, taking the Umit — > of (2 . 1 2) with Xn = nh fixed, one obtains (5.3) and the first of (5.2). 
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t .(".') ( (n.O. (1 

ol-i-i^ —•— ► oLi-^i^ . . . ► oLi-i^ — 

0,1,2, ... n~^,n,n+^ n' 0,1,2, ... n-^,n,n^-^ n' 0,1,2, ... n-^,n,n+^ 

Figure 1. The distinguished points for the eigenfunctions (^J,'', ' and (^J,^'. 



2.3. IBVPfor DLS via spectral analysis of the Lax pair 

We now use spectral methods to solve the IB VP for the DLS, namely (1.1) for n eN and f eR""", 
with q„{Qi) and qoit) given, where as before we assume q„{Qi) e {^(N) and qoit) e C(IR.q ). Before 
we do so, however, we address the issue of the well-posedness of the linear system (2.6). 

In the continuum limit, the f-part of the Lax pair evaluated at jc = depends on ^'(0, t) 
and ^,:(0, t), only one of which is given. Use of the global relation allows one to obtain the 
unknown BC in terms of the given one. In the discrete case, evaluation of the f-part of the Lax 
pair for n = requires the knowledge of q-i{t). Thus, the role of the unknown boundary datum 
in the discrete case is played by the fictitious function q-i{t). In analogy with the continuum 
limit, the solution method proceeds as though this function is given; a posteriori we will then 
show that this unknown boundary datum is determined in terms of known initial-boundary 
data via the global relation. 

A similar problem arises with Fourier methods, where one must define an appropriate 
transform so that the unknown boundary data do not appear in the expression for the solution. 
A similar situation also occurs in IB VPs for Burgers' equation [13, 14], where the solution 
depends on an unknown function that must be determined a posteriori. There, similarly to 
nonhnear PDEs solvable by the 1ST, the IB VP is reduced to a nonlinear integro-differential 
equation [14], which can be linearized for special kinds of BCs [13]. 

Eigenfunctions and analyticity. As in the continuum case [24, 26, 28], to solve the IB VP we 
consider simultaneous solutions of both the jc-part and the f-part of the Lax pair. To do this 
we again use ^„{z,t), defined in (2.7). Integrating (2.8), we then define three eigenfunctions 
uniquely determined in terms of their normalizations: namely, (p^J\z,t) for j - 1,2,3, so that 
^n\z,t) = respectively at («,f) = (0,0), as («,f) (oo,f) and at («,f) = (0, T) (cf. Fig. 1): 

4>n\z,t) = E^m(Oz""'""' +iz"/e"''"*''*'"'''(?o(f')-?-i(f')/z)df', (2.13a) 

m=0 

CO 

^?iz,t)^-Y. qm(t)z"-'"-\ (2.13b) 

m=n 

= 'e^m(Oz""'""' -/z"/e-'"<^'><'-'''(?o(f')-?-i(f')/z)df'. (2.13c) 

We introduce the domains D+ - {z e C : lm.CL>{z) < 0), which will also be convenient to 
decompose as D+ - D+^n U Z)+out, where D+^n and D+out are respectively the portions of D+ 
inside and outside the unit disk (cf. Fig. 2), namely 

D+in = {zeC: |z|< 1 Almz>0), D+om = k e C : k| > 1 Almz<0), 
D_in = {zeC: |z|< 1 Almz<0}, D_out = k e C : k| > 1 Almz>0). 

We then note that: 



IBVPsfor discrete evolution equations: DLS and IDNLS 



6 




Figure 2. (Left) The regions D+ (shaded) and D_ (white) of the z-plane where Im[tj(z)] ^ 0. 
(Right) The contours Ci,2, €2,3 and C3.1 that define the Riemann-Hilbert problem in the linear 
case (see text for details). 



(2) 

• (f>„ coincides with the eigenfunction in the IVP, hence it is analytic for |z| > 1 and 

(2) 

continuous and bounded for |z| > 1, and 0„ (z, f) ^ as z — > 00; 

• (f>lP and ^J,"^' are analytic in the punctured complex z-plane C'®^; 

• for all r > it is e''^^^" ^ as z ^ 0, 00 in D+ and e'^^^^^' ^ as z 0, 00 in D_; as a 
result, ^J,'^ and ^J,"'' ai'e bounded respectively for z e D-m and z e D+m- 

Note that (2.13) do not define <;*q '(z,f) and 4>^^^{z,i) at z = 0. In Appendix D, however, we 
compute the asymptotics of these eigenfunctions as z — > 0, and we show that (f>^^\z,t) = 0(1) 
as z ^ with Imz < and 0^ '(z, f) = 0(1) as z -> with Imz > 0. 

Jump conditions and Riemann-Hilbert problem. The difference between eigenfunctions at 
|z| = 1 and z e [-1,1] yields a scalar RHP whose solution will enable us to reconstruct the 
potential in terms of the scattering data. As before, the difference between any eigenfunctions 
solves the homogeneous version of (2.6). Evaluating these differences at («, f) = (0,0) we then 
obtain the jumps as (of course any two of the jumps uniquely determine the third one): 

0l'Vz,O-#(z,f)=z""'e-'"(*§(z,O) |z| = lAlmz<0, (2.14a) 

^*™(z,f)-<^i^'(z,0=z""'e-'"(*F(z,r) Imz = OA|z|<l, {2.Ub) 

<p^;^\z,t)-4>n\z,t)=z"-^e-"^^''>'iqiz,Q)-Fiz,T)), |z| = 1 A lmz>0, (2.14c) 

with F(z,t) = i(zfo(z,t) - f-i(z,t)), and where qiz,t) and fnizj) are respectively the z- 
transforms of the initial and boundary data; namely: 

00 ' 
§(z,f)= 2 MO/z'", /,(z,f) = /e'"<^^''Mr')dr'. (2.15) 

in=0 Q 

Note that q(z, t) is analytic for |z| > 1 and continuous and bounded for |z| > 1 , while the /„ (z, f) 
are analytic Vz and continuous and bounded for z e D+. Moreover, q(z, t) —> qoit) as z — > oo, 
while /„(z, — » as z — > 0, oo in D+. Finally, integration by parts shows that 

/„(z,f) = !z(e'"(*^„(f)-?„(0)) + (9(z2) (2.16) 

as z ^ in dD+ (i.e., along the real z-axis). As shown in Appendix B, (2. 14) are inverted by 

qn(t)^:^. § z"-'q(z,t)dz, qn(t)^^ / w'(z)e— «'/„(z,r)dz, 

\z\=l SD«„t 

(2.17) 
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for all < ? < r, where io'{z) = doj/dz and dD+out is oriented so that Rez is decreasing. 

Note that F(z,T)/z remains bounded as z ^ along the real z-axis [cf. Appendix D]. 
Thus, the RHS of (2. 14b) with n = Q does not have a pole at z = 0. The solution of the RHP 
defined by (2.14) is therefore simply obtained using standard Cauchy projectors over the unit 
circle: 

^ (2.18) 

where |^| = 1 is taken counterclockwise and dD+ is oriented so as to leave the domain to its 
left, as usual. Inserting (2.18) into the first of (2.6) then yields the reconstruction formula: 

5„(0 = -L ^ 2«-ie-''-fe)'§(z,0)dz--^^ / z"-ie-'^(^>'F(z,T)dz. (2.19) 

Of course the right-hand side of (2.19) still depends on the undetermined value q-\{t) via its 
transform /-i (z, T). We next show how to eUminate this unknown using the global relation. 

Global relation and symmetries. The global relation, which couples all initial and boundary 
values, is obtained in a similar way as in the continuum problem by integrating (2.8) around 
the edges of the domain No x [0, T] , namely for (n, ?) from (0, 0) to (0,7), from there to (00,^), 
and then to (00, 0) and back to (0, 0): 

^ 00 00 

i^e'^^\q^{t')-q-x{t')lzW ^e^^^'^' Z qm{t)lf'^'= Z ^m(0)/z"'^' . 

" (2.20) 



Equation (2.20) holds where all of its terms are defined, that is, for all |z| > 1. In terms of the 
z- transforms: 

j[2/o(z,0 -/-i(z,0] + e'"(^^'§(z,0 = §(z,0) . (2.21) 

Now note that u(z) is invariant under the transformation z — > 1 /z, and therefore so are the 
functions fn(z,t)- Moreover, z e f+out implies 1/z £ i)+in and viceversa. Hence, (2.21) with 
z ^ 1 /z gives, for all < |z| < 1 : 

i[(l/z)Mz,t)-f-i(z,t)] + e-'^^'^d/z,/) = §(l/z,0) . (2.22) 

We can then solve for f-iiz,t), obtaining, for all < |z| < 1: 

f-iiz,t) = Mz,t)/z-i(e'^'^'qil/z,t)-qil/z,0)). (2.23) 

Solution of the IBVP. Of course the RHS of (2.23) contains e''^(^'^§(l/z,r), which is (apart 
from the changes t T and z ^ 1/z) the transform of the solution we are trying to recover. 
When this terms is inserted in (2.19), however, the resulting integrand is z"~'e"^*^^^*^^~''§(l/z,/)> 
which is analytic and bounded in D+in, and whose integral over 5D+in is therefore zero. [This 
is analogous to what happens in the continuum limit; cf. section 5.] Importantly, the result 
also holds for n = 0, since decays exponentially for all ? < T as z — > in D+in. We 

then have 

9n(t) = ^. ^ z"-'e-'<^«'$(z,0)dz + ^ / z"-'^-'''^'^'[iqil/z,0)-iz-l/z)Mz,T)]dz. 

lzl=i ^ 5D+i„ ^2.24) 

Equation (2.24) provides the solution of the IBVP in Ehrenpreis form [18, 37, 32], since the 
only dependence of the RHS on n and t is via the terms ^"e"'"^^^', as in the IVP. Performing the 
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change of variable z' = l/z we can write the second term in the RHS of (2.24) as an integral 
over dD+oat- Then, since the resulting integrand, e"'"'^'^'§(z,0)/z""'"' is analytic on -D-out> for 
that portion we can deform the contour 5Z)+out onto the circle |z| = 1 and combine the result 
with the first integral in (2.24), obtaining the equivalent representation 



qnit) = :^. ^ (z"-z-")/ze-"^«'^(z,0)dz-:i / iz-l/z)z-''-'c-''"^''' Mz,T)dz, (225) 

\Z\=1 dD^onX 

where, as before, 5Z)+out is oriented so that Rez is decreasing. 



Continuum limit. The representation (2.25) is the discrete analogue of the solution in the 
continuum hmit. To see this, one can reinstate the lattice spacing h and follow the same steps 
as above. When expressed in terms of k, the solution of the IB VP then becomes: 

2 y^lh -j njh 

qn(f) = - r e-^''^'sm{nkh)q^''\k,Q)dk+- f q-^''^' sm{nkh)g{k,t)dk, 

^ ^ 

(2.26) 

where a){k) = 2(1 - cos{kh))/h^, and with 

q'-''\k, t)=hj: sminkh)q„it) , g{k, t) = 2i ^^^^ f e'^^^^'' qoit') df . 
n=i h 

It is then trivial to show that, in the hmit /j — > 0, (2.26) yield the solution of the continuum 
problem, namely (5.5). 

Remarks. Assuming existence, one can now verify that the RHS of (2.24) and (2.25) indeed 
satisfies the DDE as well as the initial and BCs. That the function defined by (2.24) solves 
the DLS equation is a trivial consequence of the fact that it is in Ehrenpreis form. When f = 
the term proportional to z~" in the first integral of (2.25) gives zero contribution, since the 
corresponding integrand is analytic, bounded for \z\ > 1, and 0(1/2""""') as z —> oo. Similarly, 
the second integral vanishes for the same reasons. The only piece left coincides with the 
RHS of the first of (2.17) at f = 0, which therefore yields the initial datum ^„(0). Finally, for 
n-Q the first integral in (2.25) is obviously zero, while the second becomes just the inversion 
integral in (2.17). Hence its result is simply qo{t). 

Even though /o(z, T) depends on values of the EC ^o(f) at all times t from to T, in 
practice (2.25) preserves causality, and the solution of the IB VP at time t does not depend on 
future values of the BCs, because one can replace T with t in (2.25). The reason is that the 
difference between the two terms is 

^ / (z-l/z)z-"-'}c-'''^'^^'-'">qo(t')dt'dz, 
dD+out t 

and Vn ^ the integrand is analytic and bounded in Z)+out. and vanishes as z ^ oo in Z)+. 
Hence, the integral is zero Vn > 0. 

For all n + Q, the second integrand in (2.25) is analytic and bounded in D-oui- Hence we 
can deform the integration contour from dD+om to \z\- 1, and substitute z ^ 1/z in half of the 
integral. The resulting expression for the solution coincides with the solution of the IB VP via 
sine series, namely (2.4). We reiterate however that (2.25) also holds for n = 0, unlike (2.4). 

Unhke sine/cosine transforms, the present method works equally well for more general 
BCs, as we show below. Also, unlike sine/cosine transforms, the present method can solve 
IB VPs for arbitrary Unear discrete evolution equations. Finally, the method can be generaUzed 
to solve IB VPs for integrable nonlinear DDEs, as we show in section 3. 
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Other boundary conditions. We now consider a IB VP for the DLS equation (1.1) in which 
the BCs are a linear combination of qoit) and q-iit) with constant coefficients, namely, when 

q-l(t)-aqQ(t)^h(t) (2.27) 

is given, a e C is a nonzero but otherwise arbitrary constant, and where in this case the labeling 
of the lattice should be such that n - -\, not n = 0, is the first lattice site. Such BCs are 
the discrete analogue of Robin-type BCs in IBVPs for PDEs, and cannot be solved using 
sine/cosine series. The present method however works equally well; the only difference from 
the previous case being that one needs to solve the global relation for a different unknown. 
Indeed, in Appendix C we show that the solution of this IB VP is given by 

qnif)^^. § z"-'e-''"fe)'^(z,0)dz--^ / z"-^e-''^«'^^dz-yaa-"-ie-'^(«>'G(l/a,f), 

(2.28) 

where 

Giz, t) = Hz - 1 /z)kz, t) + iz- a)qil/z, 0) , (2.29) 

and where Va = 1 if o- e D+out, v„ = 1 /2 if a € dD+out and Vq. = otherwise, and where the 
integral along 5Z)+in is to be taken in the principal value sense when a e dD+out- As before, 
one can easily verify that the expression in (2.28) indeed solves (1.1) and satisfies the initial 
condition and the BC (2.27). Moreover, one can also verify that, in the limit a oo with 
hit)/a = h'{t) finite, the solution of the IB VP with "Dirichlet-type" BCs [namely (2.24)], is 
recovered. 



3. Integrable discrete nonlinear Schrodinger equation 

We now turn our attention to IVBPs for the IDNLS equation (1.2). As before, we first review 
the IVP, which serves to introduce some of tools that will be used for the IB VP. We require 
the same regularity conditions on the initial-boundary data as in the linear case. 



3.1. The Ablowitz-Ladik system on the integers 

Consider the AL system (A. 4) with n e Z and f e R^, and with q„(0) given. A Lax pair for (A.4) 
is given by (A.5), where now we take <l>n(z, f) to be a 2 x 2 matrix, Q„(/) and H„(z, /) are defined 
in (A.6), and coiz) = Widnis(z) = Wdis(z^)/2, where Wdis(z) was defined in (2.3). As in the Unear 
case, we assume q„(0) e {\Z). (As in the continuum limit, the 1ST with non- vanishing BCs 
at infinity is significantly more involved, see Refs. [1, 43].) 



Jost solutions. As customary, we remove the n-dependence of the eigenfunctions as n ^ ±oo 
by introducing a modified eigenfunction as 

(t„(z,f)=yL/„(z,f)Z"e-''^(*'^3. (3.1) 

(This definition differs from the usual one by the factor e"'^^^^''^^^ which has been added for 
consistency with the the IB VP, discussed in section 3.2. With this choice, the scattering matrix 
will be independent of time.) Then /z„(z, t) satisfies the following modified Lax pair: 

yU„+l-ZyU„=Q„jU„Z"', fln + i<^(z)[o-3,ldn]^y^nl^n, (3.2) 

where ZA = ZAZ"' It is also useful to use the integrating factor e'^^^iA) = Q^Sa-i j^^-iea-^ 
(cf. Appendix A). Then, the function 

^n{z,t) = Z-"e'^'^'^^ix„iz,t) , (3.3) 
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solves 

T„+i-1'„ = Z-iZ-V"(^)"^3(Q^)<I/^. *„=Z-V'^(^>"^3(H„)>P„. (3.4) 

One can now easily "integrate" (3.4) and thereby obtain the solutions of (3.2) which reduce 
to the identity matrix as n — > +oo: 

tl-;^\z,t) = I + "2 Z"-'"(Q^/iL'') , ti\z,t) = I - Z-i 2 Z"-'"(Qm/iL") • (3.5) 

m=—oo m=n 

Of course, unUke the linear case the eigenfunctions are now defined in terms of summation 
equations (the discrete analogue of integral equations). 

As in the linear problem, (3.5) imply certain analyticity properties for the eigenfunctions. 
More precisely, let jj.n\z,t) = (j/J'^Ki/J''^^), j - 1,2, where the colunm vectors Hn'^\z,t) and 
l4i''^\z,t) denote respectively the first and second column of i4Piz,t). These columns are 
analytic in the following regions [6]: 

\z\>l, IzKl, 

Moreover, these columns are continuous and bounded on the closure of these domains. 
These properties innmediately yield those of o'-J^izJ) = ju^-'\z,0Z"e"''^(^'"^3 for j = 1,2: 
^^n'^\z,t) and ^^n'^\z,t) are analytic for |z| > 1, and ^^!'^\z,t) and ^^n'^\z,t) for |z| < 1. 

Scattering matrix. Equation (A.Sa) implies det (S>„+\ = (1 - q„p„)det <&„. Therefore 

det <^ = n (1 - qmPm) , det 0)1'' = n (1 - qmPm)-' =: 1/C„ . (3.6) 

m~—oo m—n 

(Note detO„ - det//„.) Equations (3.6) mark a significant difference of the discrete case from 
the continuum case, where such determinants are independent of both the potential and the 
independent variable (cf. section 5). 

For the focusing IDNLS [namely, (A.4) with p„ - vq^ and v = -1], 1 -q„p„ - 1+ \qn\^, 
and therefore detj/J^ Vn e Z for j = 1,2. For the defocusing case (v = 1), however, it is 
necessary to assume that \q„\ 1 Vn € Z in order that det/ij/' to be guaranteed to be nonzero. 
Hereafter we will assume that qnPn ^ 1 Vn e Z. Moreover, we will require that the product 

OO 

C-oc = det 0)^' = 1 / det (D?i =0(1- qnPn) 

n=—oo 

be finite, which will simplify the study of the scattering coefficients. Under these hypotheses, 
the matrices Oj,'^ and 3>„ are both fundamental solutions of the scattering problem (A.5a). 
Hence they must be proportional to each other: <S>^^\z,t) = <S>^^\z,t)A(z) on |z| = 1, where 
A(z) - {ajfiz)) is the 2x2 scattering matrix. In terms of the modified eigenfunctions: 

fii'\z,t) = £\z,t)Z" c-'^^'^'^^Mz) . (3.7) 
Or, in component form, 

/ii''^'(z,f) = an(z)fi^n'''\z,t) + z-^''e^^'^'a2i(z)fif''\z,t), (3.8a) 
li-^hzj) = z'"^-^''^'^'an(z)n'^''^\z,t) + a22iz)nf'''\z,t). (3.8fc) 
The above relations imply A(z) = lim„^ccZ-«e'"(**3^ii)(2,;) = hm„^^>p^i'(z,?), that is, 

A(z) = l + Z-i 2 Z-"s-''^'^'^^iQnit)n';^\z,t)). (3.9) 



IBVPsfor discrete evolution equations: DLS and IDNLS 



11 



The scattering matrix A(z) is independent of time, since A(z) = Um„^ooYi^\z,0, and 
lim„^oo*i'n(z.O = 0. Equation (3.7) also implies det A(z) = det <l)^'(z,f) = C_oo, as well as 

The analyticity of the eigenfunctions then implies that a\i{z) and a^iiz) can be analytically 
continued off the unit circle, respectively into the domains |z| > 1 and \z\ < 1, but ai2(z) and 
aiiiz) cannot. It is also useful to introduce the reflection coefficients 

Pi(z) = «2i(z)/«ii(z), P2iz) = auiz)/a22iz). (3.11) 



Symmetries. When pn(t) - vq*(t), the scattering problem (A. 5a) admits an important 
involution, which can be conveniently written introducing the matrix o-y defined in (A.lla). 
Indeed, when p„ = vq*, if <l)„(z,f) is a solution of (A.5a), so is the matrix 

<D;(z,f) = (^v<l>:(l/z*,0, (3.12) 

Then, comparing the asymptotic behavior of the first and second columns of the Jost 
eigenfunctions as n — > ±oo one obtains, for 7 = 1,2, 

a>i-''^'(z,0 = crv(a)i^'*>(l/z*,0)*, o';^'''\z,t) = y(r,(a>2'^'(l/z*,0)*. (3.13) 
The above relations imply the foUowing symmetries for the elements of the scattering matrix: 

a22(z) = an(l/z*), a2iiz) = va\2il/z*). (3.14) 
In turn, these imply piiz) = vpj(l /z*). 



Discrete spectrum. The proper eigenvalues of the scattering problem (A.5a) are the values 
z = Zj with \zj\ < 1 and z = zj with > 1 for which there exist eigenfunctions bounded 
Vn e Z. From the asymptotic behavior of the Jost solutions one can see that such eigenvalues 
occur whenever the appropriate left- and right-sided Jost solutions are proportional, namely 
^^l-^\zj,t) = bf^^f-^izj,t) and ^^n''^\zj,t) = bf^f-^\zj,t), or equivalently: 

/^i''"\z;sO = 5fzf''e2'--C^Vf«\z,,0, £''^(zj,t)^bfz^^'^-^^^'fe'%j,t). (3.15) 

The Wronskian representations (3.10) then imply that such eigenvalues are the zeros of the 
scattering coefficients: aii(z/) - and fl22(Z;) = 0, respectively. (As in Ref. [6] we assume 
that ajjiz) + for all |z| = 1 .) Since no accumulation points of such zeros can exist (because of 
the sectional analyticity of the scattering coefficients), it follows that there is a finite number 
of them. As in Ref. [6] we assume all of these zeros are simple. (The case of multiple zeros 
can be studied as the coalescence of simple zeros, by analogy with the continuum case [44].) 
Since ajj{z) are even functions [6], z = zj is a zero of a22(z) iS z = -zj is, and similarly for 
aii(z). Moreover, the symmetries (3.14) imply that z = Zj is a zero of 022(2) ifl" Zj = l/z*j is a 
zero of aii(z). Thus, discrete eigenvalues appear in quartets. 

The inverse problem wiU involve the modified eigenfunctions i/^'^\z,t)/aniz) and 
Hn''^\z,t)/a22iz). Equations (3.15) imply 



Res 



bfzf"c''''^'jVn'^\zj,t), Res^^ii-^ ^bjzfc-'^'^Vn''\zj,t), 
■> 7=7.- a22(Z) J 

(3.16) 



z=zj L aii(z) 

where bj = bfla'^^izj) and bj = bf^/a[^(zj) are referred to as the norming constants. The 
symmetries of the scattering problem imply bj = -vibj/z^)*. 
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Asymptotics. The asymptotic behavior of the eigenfunctions as z — > or z — > cxj can be 
obtained from (3.5). For example, for iJ'^\z, t) it is 

lIn\z,t) = \+Qn-i'Zr^ +0{Z-^) asz^(co,0), (3.17) 

where z — > (zl,zr) indicates z — > zz, in the first column and z ^ Zr in the second one, and the 
asymptotics corresponding to 0(Z'") is defined in Appendix A. Equation (3.17) will allow us 
to reconstruct the potentials from the asymptotic behavior of f/^^: 

Q„it)= lim (//(i^\(z,0-l)Z. 

Z->(oo,0) "^^ 

The asymptotic behavior of iJn\z,t) is obtained in a slightly different way as that of i/n\z,t), 
and the result is also different. More precisely, in Appendix D we show that 

Cnf/n\z,t)^\-QnZ + 0(Z^) as Z ^ (0, oo) . (3.18) 

Also, inserting (3.17) into the diagonal elements of (3.9) one obtains the asymptotic behavior 
of the analytic scattering coefficients: 

1 

an(z)^l + -j Z q„(t)pnit) + Oil/z), asz^oo. 

Z «=-co 

which by synometry also determines the behavior of a22(z) as z — > 0. 

Inverse problem. The inverse problem is the RHP defined by (3.8) for |z| = 1: 

-£''\z,t) = z-^«e2'-^«'pi(z,0A.i^^^\z,0, (3.19a) 

aniz) 

^^^^-/.i2'«>(z,0=z^"e-2'-(*p2(z,0^^^^^^ (3.19fc) 

(2 L) 

where pi(z) and piiz) as in (3.11). Unlike the continuum case, the asymptotics of //„ ' (z,?) 
as z — > 00 depends on the values of the potentials qmit) and Pmit) for all m > n through C„ 
[cf. (3.6)] . This problem can be circumvented by introducing the following renormalizations: 

,(1 



J«(z,r) = 



The matrices MJ(z,f) are sectionaUy meromorphic for |z| < 1 and |z| > 1, respectively. 
Moreover, (3.19) yields the following jump condition for the matrices M^(z,0 on |z| = 1: 

m;(z,/) = m;^(z,/)(i-J«(z,0), (3.20) 

where the jump matrix J„(z, t) is 

Piiz)P2(z) z2"e-2-b)'p2(z) ■ 
-z-2"e2'"W'pi(z) 

Moreover, M^(z,f) have the following asymptotic behavior: 

M-(z,0 = l + -fn ~'^"n'^"] + 0(l/z^) asz^oo, (3.21a) 

^«^^'^> = ('o''" C„yi-l 'VyOiz') asz-O. (3.21.) 

In the absence of a discrete spectrum [that is, if aii(z,0 ^ for |z| > 1 and a22 for |z| < 1] 
the matrix functions M^ix,t,k) are analytic in their respective domains. 
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In particular, (3.21a) allows the RHP (3.20) to be solved via the Cauchy projectors 
over the unit circle, as in the linear case. Of course, unlike the linear case the solution is now 
expressed in terms of a matrix integral equation: 

Kiz,t) = \ + ^. / M:(^,0^^^d^. (3.22) 

If 1=1 b Z 

The asymptotic behavior of \\^^{z,t) as z -> is easily obtained from (3.22): 

m:(z,/) = I + :^ / m:(^,OJ«(^,0^ + :^ / Kii,t)Ui,t)% + 0{z\ 

(3.23) 



2^' If 1=1 i ^^^\(\=i ^' 



Comparing the Umit as z ^ of (3.23) with (3.21fo), we see that the off-diagonal portion of 
the first integral in (3.23) is zero, a fact which is not entirely obvious otherwise. (This integral 
is missing in the corresponding formula in Ref. [6].) Then, comparing the (1,2) components 
of (3.23) and (3.21fc) we obtain the reconstruction formula for the solution of the IVP: 

qn(t) = :^. J z'"e-2-feV2(z)(/z®i(z,f))ndz. 
^'^^ lzl=l 

Linear limit As in the continuum limit, the 1ST is the nonlinear analogue of the linear 
transform pair. Namely, if Q„ = 0(e), then /z^^^ = I + 0(e) and 

oo 

A(z) = l+Z-i 2 Z-"e-''^(^^"^3Q„(0 + O(e2). 

n=— OO 

Thus 

CO 1 

PliO = 2 ^-'"-ie-2'"(f>'^„(0 + 0(e2) = -§(^^ 0) + 0(e2) , 

n=-oo ^ 

where q{z, t) is the Unear z-transform defined in (2. 1). Similarly, 

qn{t)=^. J ^'"e-2''"(«'p2(0d^ + O(e2) = -L J z"-'c-''^i^^(^'>'q(z,0)dz + O(s^), (3.24) 
|f|=l |z|=l 

where the change of variable = z was performed in the RHS of (3.24), and where 
^'>idnis(^) = j'^disC^^)) as discussed in section 2.1. 

3.2. The Ablowitz-Ladik system on the naturals 

We now consider the IB VP for the IDNLS. That is, we solve (1 .2) with n e N, ? e R+ and with 
q„(0) and qo(t) given. The approach we will follow is a combination of the method for the 
IVP for the IDNLS on the integers and that for the IB VP for the DLS on the naturals. 

Eigenfunctions and analyticity. Making use of the modified eigenfunction ^^^(z,?) in (3.3), 
we define three eigenfunctions Hn\z,t) which reduce to the identity matrix respectively when 
{n,t) = (0,0), as (n,f) (oo,t) and at (n,t) = (0,T): 

/iii\z,0 = l+Z-i"2 Z"--(Q^(0/iil'(z,0) + z4e""^'^^'"''^^^ (3.25a) 

m=0 

£\z,t) = I -Z-i 2 Z"-"'(Q^(t)t/^\z,t)), {3.25b) 

m=n 

fi'^\z, = I + Z-1 "2 Z''-'"iQ^it)f/^\z, t)) - Z" Je--fe)(^-'')-3 (Ho(2, t')fi'^\z, t')) dt' . (3.25c) 

m=0 t 
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Figure 3. (Left) The regions D+ (shaded) and D_ (white) of the z-plane where Im[(x)(z)] ^ 
in the nonhnear case, with D± = D+i^ U D+om ■ (Right) The contours Li , . . . , L4 that define the 
Riemann-Hilbert problem (see text for details) 



Note that /j„ (z,t) coincides with the eigenfunction in the IVP, defined in (3.5). As in the 
linear case, we partition the complex z-plane into the domains D+ defined as D+ = {z e C : 
Ima;(z) < 0}. We then write D+ - £>+in UD+out where the subscripts "in" and "out" denote the 
portions of D± inside and outside the unit disk, respectively. That is (cf. Fig 3), 

D+in = {z e C : |z| < 1 A arg z € (0, n/l) U (tt, 37r/2) ) , 
D_in = {z e C : |z| < 1 A ai-gz 6 (7r/2,;r) U (3;r/2,27r)}, 
D+out = {z e C : |z| > 1 A argz e (7t/2,7t) U (37r/2,27r)) , 
D_out = {z e C : |z| > 1 A argz e (0,n/2) U (tt, 37r/2)) . 

Then, in a similar way as in the IVP on the whole line and the IB VP in the linear problem, we 
can obtain the regions of analyticity and boundedness of the eigenfunctions. More precisely, 
writing again f/J\z,t) - if/J'^^f/J''^^), we have: 

• fi,i\z, t) and fifi\z, are analytic in the punctured complex z-plane C ^'^^ ; 

• il^'^\z,t) is continuous and bounded in D+oui, 

• the restriction of jXn'^HzJ) to Z)_in is continuous and bounded in D-{^; 

• yujj ' '(z,f) is continuous and bounded in Z)_out; 

• the restriction of ji), ' {z,t) to D+m is continuous and bounded in D+m; 

• ^„ ' (z,f) is analytic for |z| < 1 and continuous and bounded for |z| < 1; 

• yu[f'^\z,/) is analytic for |z| > 1 and continuous and bounded for |z| > 1. 

The analyticity of the eigenfunctions is formally proven via Neumann series as in the IVP [6] 
and as in the IVP for the continuum case [29]. However, showing the continuity of LL'^'^\z,t) 
and ii„ ' izJ) at z = is nontrivial, and it requires studying the asymptotic behavior of the 
eigenfunctions as z — » (see Appendix D). 

Scattering matrices. The relation det <I>„+i = (1 - qnPn) det 0„ still holds. Therefore det <b^p 
and detO*f^ are still given by (3.6), and detO^'' - detO*f\ [Note that ji, - L^i+^VK 
impUes (det//), = tr(L + R) det/i, and in our case both L and R are traceless; cf. (A. 5) 
and Appendix A.] Hence, under the same regularity hypotheses as before, O^^^ and <^^^^ 
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are each fundamental solutions of the Lax pair (A.5). We can therefore write the following 
relations among the modified eigenfunctions: 

l^?iz,t) =/i<'\z,/)Z"e-'"<*'^3s(z), (3.27a) 

/if = /<i'>(z,0Z"e-'^<^)"^3S(z,r), (3.27/7) 

which hold wherever all terms are defined, namely: the first column of (3.27a) holds for 
< Izl < 1, the second column for |z| > 1 and {3.21b) holds Vz 0. Thus 

s(z) = il^\z,0) , S(z, T) = (e'"«^*3^[,'\z, T))'^ . (3.28) 
Equation (3.28) allows us to write integral representations for the scattering matrices: 

s(z) = I - Z-i 2 Z-"(Q„(0/*P(z,0)) , (3.29a) 

n=0 

r 

S-\zJ) = I + / e'^'^'^^{Hf,{z,t)nf{z,t))At. (3.29fo) 



Note that s(z) is again independent of time, since s~^{z) = \imn-^ooZ~"e"^^^^'"'^iJn\z,t) = 
as in the IVP. Note also that (3.27) imphes 

det s(z) = 1 /Co , det S(z, T) = l. (3 .30) 

The analyticity properties of iJ'^\z,t) are the same as those of /J^\z,t)- However, i/^\z,t) 
enjoys larger domains of analyticity and boundedness than i/^\z,t). The analyticity and 
boundedness regions of the scattering matrices are determined correspondingly via (3.28): 

• Sl(z) is analytic for |z| < 1 and continuous and bounded for |z| < 1 ; while S/f(z) is analytic 
in |z| > 1 and continuous and bounded for |z| > 1 ; 

• S(z,T) is analytic in C^**J; moreover, SLiz,T) is continuous and bounded in D-, while 
Sk(z, T) is continuous and bounded in D+. 

The above boundedness properties of S(z, T) can be obtained as follows. Let us write the 
matrix S(z, T) as 

A(z,r) B(z,T)' 



S(z,T) 15(^7^) Aiz,T) 

As we show below, the symmetries of the problem imply that A{z, T) and B{z, T) can be 
obtained respectively in terms of A{z,T) and B(z,T). Hence, we only need to discuss the 
properties of A(z,T) and B(z,T). Recall that S{z,T) is an entire function of z, and note 
that (3.30) implies 



vz,i) y-B(z,T) A(z,T) }■ 



Then (3.28) and the analyticity properties of ij.q \z,T) imply that A(z,T) is bounded in D-. 
Also, (3.29Z>) and the integral representation (3.25a) with n = can be used to write a 
Neumann series for S~Hz, T), which in turn can be used to prove analyticity and boundedness 
ofB(z, r)inD_. 

The involution symmetry discussed when dealing with the IVP is a local property. 
Therefore, when /5„(f) = vq^{t), it also apphes for the IBVP. That is, (3.13) still holds, as 
does (5.30) for j = 1,2,3. This implies 

/a(z) vb*il/z*)\ ^.T..(Aiz,T) vB*il/z*,T)\ 
^^^'~\biz) a*il/zl j' "^^^''''[Biz,!) A*il/z*,T) j' ^^-^'^ 
Note that (3.30) imply 

a(z)a*(l/z')-vMz)^*(l/z') = 1/Co, 
A(z,7')A*(l/z*,T)-vB(z,T)B*(l/z*,T) = 1. 
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Asymptotics. Since //„ (z,t) coincides with (3.5), its asymptotics as z — > (0,oo) is still given 
by (3.18). Also, in Appendix D we show that, even though the definition of fi^\z,t) and 
l^^n^zj) involves time integrals, it is still n^\z,t) = l + 0(Z~') as z — > (oo,0) for 7 = 1 and 7 = 3 
in their respective domains of boundedness. More precisely, for all n > it is 

H^'^{z,t) = \+Qn-\{t)Z-^+0{Z-^), asz^(oo,0) (3.32a) 

for j = 1,3, and the limits are restricted the appropriate regions of the complex plane, where 
the corresponding columns are bounded. For n = it is instead 

nf{z, = 1 + (Q-i (0 - e-''^(^'"^3Q_j(o)) z-i + 0(Z-2) , {332b) 

nf{z,t) = I + (Q-i(f)-e-'"(^^('-^)*3Q_j(r))Z-i +(9(Z-2), (3.32c) 

as z — > (00, 0). The above yield, for all n > 0, 

Q„-i(0= lim {^Ii\z,t)-\)Z, for 7 = 1,3. (3.33) 

Z->(oo,0) 

Also, the asymptotic behavior of the eigenfunctions determines that of the scattering matrices. 
In particular, from the second of (3.28) we have 

A*(l/z*,r) = 1 +0(l/z^), B*{\lz*,T) = 0(l/z) as z ^ c» in 5+out, (3.34a) 
while (3. 27b) implies 

A*(l/z*,7') = l + 0(z2), B*il/z*,T) = Oiz) asz^OinD+in. (3.34fc) 

Similarly, (3.18) and (3.27a) yield 

a*(l/z*) = l/Co + 0(l/z^), fo*(l/z*) = 0(l/z) asz^c»in5+out. 

(3.34c) 



Riemann-Hilbert problem, solution and reconstruction formula. We now formulate the RHP 
whose solution will enable us to obtain a representation for the solution of the AL system on 
the naturals. For later reference, we introduce the quantities 



r(z) = 



vb*{z) 



R{z,t) = 



B*(l/z*,t) 
A*(l/z*,0 ' 



r(z) = 



B(z,T) 



a*(l/z*)d*(l/z*) ' 



with 



d(z) = a{z)A\\lz*,T)-vb{z)B\\lz*,T) . 
Note that R{z, T) is defined Vz e C except where A*(l /z*,T) = 0, r(z) is defined for z e L3 UL4, 
diz) for z 6 D±m, and yiz) for |z| = 1. Moreover, d\\/z*) = 1/Co + 0(l/z^) as z — > 00. In the 
analysis of Unearizable BCs, it will be useful to write r*(l/z*) in terms of only a(z), b{z) and 
Riz,T) as 

, ^ T) 

a(z)ia(z)-vb(z)R(z,T))' 
Finally, we introduce the normahzation matrix C„ = diag(l /Co, C„) . 

We are now ready to formulate the RHP, which we do using (3.27). We introduce the 
matrix functions M^(z,0 defined as: 



K(z,t) ■ 



M-(z,t) ■ 



c 



d(z) 



(z,t) 



(l,R) 

I 

(i) 



z e £)+in, 
z e D+out, 



/ (2,L) ^\ 

"T" 'fl(z)i' 

^"U*(l/Z*)'^" I' 



z e £)_in , 

Z G D-out . 



(3.35a) 



(3.35fc) 
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Note that M*(z,?) are sectionally meromorphic respectively for z € D+ and z e Z)_. Moreover, 
after some tedious but straightforward algebra, equations (3.27) yield the jump conditions as 

M,7(z,f) = M;^(z,f)(l-J«fcO), z€L, (3.36) 
where the contours L- L\ U L2 U L3 U L4 are (cf. Fig. 3) 

and the jump matrices J^^\ . . . , J*,"*^ are defined by 

J«(.o = (2 --e-7*(l/^^^)), 

Jlf\z,f) = I - (I -Jii>)(l -Ji'Vkl-jf ), z^U. 
As in the IVP, we first consider the case in which no discrete spectrum is present. For the 
IB VP, this corresponds to assuming that a(z) ^ for z G f-in and d{z) for z e D+in. In this 
case, the matrix functions MJ(z,f) are analytic in their respective domains. Also, M„(z,0 — > I 
as z ^ 00 thanks to (3.18), (3.32a), (3.34) and (3.35). Hence the matrix RHP (3.36) is solved 
by the Cauchy projectors over the contour L, namely = l/(27r0 /^[l/(^'-^)]dA:'. That 
is, 

M:(z,0 = l + ^/M:(^,0^^|3j^df. (3.37) 
Equation (3.37) also yields the asymptotic expansion of M;5^(z, f) as z ^ 0, namely, 

M:(z,0 = I + ;^ /M:(f,OJ„(f,f)^ + -^ /M:(f,0J„(^,f)^ + (9(z2). 

Note that we can write (3.38) as 

M;^(z,0 = diag[l/(CoC„),CoC„] + / M+(f,0(l - J„(^,0) ^ + ^(z'). 

^ ^ (3.39) 

Now note that the matrix C~'M„(z,0 satisfies the n-part of the Lax pair (A.5a). Also, thanks 
to (3.18), (3.32a) and (3.34), it is 

M+(z, t) = diag[l /C„, Co] + 0(z) as z ^ . 

Hence, substituting the asymptotic expansion of M„(z, f) into (A.5a) and comparing the (1,2)- 
components of the 0(z) terms, we can recover the scattering potentials as 

<7„(f) = lim(M,;^i(z,f)-l)i2/z- (3.40) 

Taking the (l,2)-component of (3.39) and comparing with (3.40), we then obtain the 
reconstruction formula for the solution of the IDNLS equation on the natural numbers: 

qn(t) = / z'"e-2-»y(z)(M:^i(z,0)ii dz + / z'"e-2-»r (1/z*)(M:^i(z,0)ii dz 



+ ^.JiyCoMzf-Co + l)(M„^i(z,0)i2 f + ^ (i - ^) /(M:.ifc0)i2 1 

+ / z'"e-2-«'r*(l/z*)(M:^i(z,0)iidz. (3.41) 

2m Co 
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Global relation. As in the linear problem and the continuum limit, the unknown boundary 
datum can be obtained in terms of the known initial-boundary conditions using the global 
relation and the syrmnetries of the system. 

Integrating (3.4) around the boundary of the region No x [0, i], one obtains 

Z / e''^fe>'^3(Hofcj')juofcf'))df' + e''"(^)"^3 ^ Z-«(Q„(OAi«fcO) = 2 Z-«(Q„(OKfcO)). 

Q n=0 n=0 

(3.42) 

When (3.42) is evaluated with jJ.„(z,t) = yuj, (z,t) and t - T, the first and second columns of 
the resulting equation are valid respectively for z e D+i^ and z e 5±out- Moreover, the RHS 
of (3.42) becomes Z(l -s(z)) thanks to (3.29fl). Finally, using (3.27), we can write the first term 
and the second term in (3.42) respectively as Z(S"'(z, T) - l)s(z) and Z(l -e''^(^^"^3/i[)^\z,r)). 
We therefore have the following global relation in terms of the scattering data: 

S-'(z, T)s{z) = I - e'"(^^^'^3G(z, T) , (3.43) 

where 

oo 

Giz,t) = Z-'Y.Z-"iQ„fi\z,t)). 

n=0 

Like for (3.42), the first and second colunrn of (3.43) are respectively valid for \z\ < 1 and 

(2) 

\z\> 1. Also, from the analyticity domains of (z„t) it follows that Giizj) is analytic in 
|z| < 1 and GR(z,t) is analytic in |z| > 1. Taking the (1,2) component of (3.43) we have 

A*(l/z*,T)b*(l/z*)-B*(l/z*,T)a*(l/z*) = -ve^'^^'^^ G(z,T) , \z\ > l,^^ 

where 

G{Z,T) = i Z-^"-'qn(T)(f/n\z,T))22- 
n=0 

Also note that the RHS of (3.44) is bounded for z e D+out- Then, for z e D+out the RHS 
vanishes in the limit r — > oo, implying 

A*(l/z*,T)b'(l/z*)-B*(\/z*,T)a*(l/z*) = 0, zefi+out- (3.45) 

For finite values of T, letting r(z) = b(z)/a(z), the global relation is now 

B*(l/z*,r)-r(l/z*)A*(l/z*,r) = ve2'"«^G(z,r)/fl*(l/z*). 

Since Giz,t) = 0(l/z) as z — > for z e D+oau multiplying by zc-^"^(dt integrating over 
^^+out [where D+out = {O+out Almz > 0}] we obtain the integral relation: 

/ ze-^''^''>'(B*il/z*,T)-ril/z*)A*il/z*,T))6z = 0. (3.46) 

5D+OU1 

This is the discrete analogue of the one that in the continuum case is used to obtain the 
Dirichlet-to-Neumann map [28]. In the discrete case, however, the unknown boundary datum 
can be obtained using an alternative, simpler method, as we will show in section 4. 



Linear limit. The linear limit of the solution (3.41) of the IB VP for the IDNLS equation 
coincides with the solution of the IBVP for the DLS equation, as we show next. Suppose 
Q„(;) = 0(e). From (3.25) it follows that yu„ = I + 0(e). RecalUng (3.29a), we obtain, to 0(e): 

y{z) = ~q{^,Q), d{z) = l, Co = l, r{\lz*) = iv{^J-i{^,T)-zm^,T)]. 
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Thus (3.41) yields, to 0(e), 

qnit)=^. J ^2«-ie-'"(f'>'^(f^0)d^+;i / ^2«e— (f')'(i/_i(^2,r)-^/o(^^r))d^. 

If 1=1 L^+U ' 

where the integrals are taken in Cauchy's principal value sense. Now note that the contour 
Li UL4 can be deformed to c?Z)+in by Cauchy's theorem. Performing the change of variable 
= z, we then obtain that the linear Umit of (3.41) coincides with the solution of the IB VP 
for the DLS on the natural numbers, namely (2.24). 



Continuum limit. Reinstating the lattice spacing h, it is easy to show that the Lax pair for 
the NLS is the continuum limit of that for the AL system as h [4, 6]. The continuum 
limit is formally obtained by writing the solution of the discrete case as Qn(t) = hq(nh,t) and 
Pnit) = hp{nh,t). Then for z - s'^'', the Lax pair (3.2) becomes 

- ik[cr3,ti„] = Q„(f)M„ + 0{h^) , 

h 

fin + ia>(k)[o-3,pn] = H„(f,^)jU„ + 0(h) , 

where now (oik) = (1 -cos2kh)/h^, with p„ = p{nh,t,k), qn = q{nh,t) and pn = p{nh,t) for 
brevity, and where 

-qnPn-i iqn-qn-i)/h + ikiq„ + q„-i) 

-iPn- Pn-\)lh + ik{pn + Pn-\) qn-iPn 

Correspondingly, the Jost solutions are obtained from (3.25), for example, 

n';l\k,t) = I +/j "2 e*''("-'">*3(Qm(0/^m (A;,?)) + je'^''''''-'^''^^'-''^^^HHiO,k,t')p'^\k,t'))dt' . 

m=0 Q 

As ft ^ with X = nh fixed, we have a)(k) 2k^, together with H„(f,fc) — > H(x,t,k) and 
l4Pik,t) ilj\x,t,k), j = 1,2,3, where p^j\x,t,k) are the Jost solutions for the IB VP of the 
NLS, namely (5.38). Note also that C„ — > 1 as /i — > 0. Hence, in the continuum hmit, the 
solution of the IB VP for the IDNLS becomes exactly that of the IB VP for NLS. 

The result can also be verified directly via the continuum hmit of the solution (3.41). 
Explicitly, since Co = 1 + 0{h}), as ft —> we have 



Hn{t,k)^i 



Qn{t) = -^. / ^2"e-2-(f)V(^)(M:,i(f,0)ndf 



2m 1^1^ J 

+ ^ / ^2"e-2-(f)'r(l/r)(M;,i(^,0)„d^+ ^ / lr(^)|^(M;,i(^,0)i2 ^ + 0(ft2). 

The oriented contour L\ U L4 can be deformed onto |^| = 1 since the corresponding integrand 
is analytic in D-m. In terms of q{nh, t) = Qn{t)lh, and performing the substitution ^ = e''^'', we 
then have 

q{nh,t) = -- / e2*'*^''-'^(*»r('fc)(M^+i(A;,0)iie''=''dA; 

^ -nIh 

n/h h "^'^ 



- J ye2'("*''-'"»)r(r)(M;J^i(^,f))iie'*'M^-- / v|r(^)|2(M;^i(^,/))i2e-*'dA: 

^ n/2h 

h "'^'^ 

- / y\7(-k)\\M;^,i-k,t)\2c'"'dk. (3.47) 



^ -n/h ^ nl2h 

h "'^^ 
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Now note that, since Q~'^'-'^'T*{k*){U^^^{x,k,t))u is analytic and bounded for (Rsk & 
[-n I h,Q]) f\(}mk> 0) [which becomes Cn in the limit/; — > 0], the portion of the corresponding 
integral on the negative real axis can be deformed onto the positive imaginary axis. Then, 
taking the continuum Umit of all the integrals in (3.47) we obtain that q{nh,t) coincides with 
the solution of the IB VP for the MLS, namely (5.45), in the Umit h^O. 



Remarks. A few comments are now in order: 

• Equation (3.41) provides the Ehrenpreis [18, 37, 32] representation for the solution of 
the IB VP for the IDNLS, in analogy with Ref . [27] in the continuum Umit. 

• One can now use (3.41) as a starting point to formally prove that the function q(x, f) given 
by the reconstruction formula satisfies (1.2) as well as the initial-boundary conditions, 
using the dressing method, as in Ref. [29] in the continuum limit. 

• In the continuum problem, the location of the jumps is the union of the jumps for the 
scattering problem in the linear case and those of its adjoint. In the discrete problem, 
however, this is not the case. Indeed, the extra jump along the imaginary axis arises as a 
consequence of the rescaling z — > when going from the linear to the nonUnear case. 

• The scattering matrix S(z, T) involves T explicitly. In Appendix E, however, we show 
that the solution of the IB VP for the AL system on the naturals does not depend on future 
values of the boundary datum. 

• With the due modifications, the method presented here can also be used to solve the 
IB VP for all members of the Ablowitz-Ladik hierarchy. Moreover, the method can be 
generaUzed to any integrable differential-difFerence evolution equation. 



4. Elimination of tlie unknown boundary datum, linearizable BCs and soliton solutions 

4.1. Elimination of the unknown boundary datum 

The scattering matrix S(z, T) depends on the both the known and the unknown boundary 
datum. In the linear problem, it was possible to overcome this difficulty by making use of the 
fact that the transformation z — > 1/z leaves the transforms of the boundary data unchanged. In 
the nonlinear problem, however, the matrix S(z, 7") is not invariant under this transformation, 
because it is defined in terms of the eigenfunction iJ^\z,t), which is not invariant under 
z—^ l/z. As in the continuum case [28], the determination of the unknown boundary datum 
in terms of the known initial-boundary conditions is in general a nontrivial issue. 

For linearizable BCs it is possible to express the RHP only in terms of the initial data, 
as we show in section 4.2. This is not possible for generic BCs, however. In this case one 
must solve a coupled system of nonUnear ordinary differential equations (ODEs) to obtain 
simultaneously the unknown boundary datum q-iit) as well as scattering coefficients A{z, T) 
and B{z, T), as we show next. 

The boundary data enters the RHP only via the ratio Riz, T) = B(z, T)/A{z, T) appearing 
in r(z). Recalhng (3.1) and (3.28), we have S(z,t) = ^-'^{z,t)t-''^''^"^\ where the matrix 

4.fz - a.(iVz - ( -v^-'^^^^^d/z*) \ 

<^{z,t)-%{z,t)-y _^Uo(z)tB(z) e'^fe)'A(z) j' ^^"^^ 
satisfies the ?-part of the Lax pair (A.5) for n = 0, namely: 

«!) = (- iajiz)o-3 + Ho(z, t)) 3) , (4.2) 
together with the initial condition 0{z, 0) = I . 
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The term Ho(z,0 in (A.6fo) contains q-\{t), of course. Note however that using (3.33) 
with n = 0, we can express q-\{t) in terms of ^z,?): 

^_ia) = lim(//<,'\z,0)i2/z zeD_in. (4.3) 

The simultaneous solution of (4. 1) and (4.3) provides the unknown boundary datum as well as 
the auxiliary spectral functions A{z, t) and B{z, t), allowing one to completely define the RHP 
and therefore we also obtain the solution of the inverse problem. Note that this procedure is 
significantly simpler than that requried to obtain the generahzed Dirichlet-to-Neumann map 
in the continuum case [28]. 



4.2. Linearizable boundary conditions 

Like in the continuum case, there is a class of BCs, called linearizable, for which it is possible 
to obtain the unknown boundary datum via only algebraic manipulations of the global relation. 

Recall that A(z,0 and B{z,t) are given in terms of 0(z,O = ju[,''e"''^'^^'"^3 by (4.1) which 
solves the ODE (4.2) together with the initial condition 0(z,O) = I. Since w(l/z) = w(z), 
the matrix 3)(l/z, f) satisfies equations identical to (4.2) except that Ho(z,0 is replaced 
by Ho(l/z,f)- If there exists a time-independent matrix N(z) such that 

N(z) ( - iu{z)(Ti + Ho(z, 0) = ( - iojiz)(r3 + Ho(l /z, t)) N(z) , (4.4) 

it is then easy to show that 

4)(l/z,0 = N(z)d)(z,0N(z)-i. (4.5) 

A necessary condition for (4.4) to be satisfied is obviously that det[-/w(z)cr3 + Ho(z,f)] = 
[(z^ - 1 /z^)iqoP-i - Q-iPo)]^ be invariant under the transformation z — > 1/z. In turn, for this 
condition to be satisfied one needs 

qoP-i-q-ipo = 0. (4.6) 

In the reduction p„(t) - vq*(t) to IDNLS, (4.6) is satisfied by the discrete analogue of 
homogeneous Robin BCs: 

q-l-Xqo^O, ;r6R- (4.7) 
These BCs had been previously identified via algebraic methods [31]. For the BCs (4.7), we 
can solve the system (4.4) for N(z), obtaining Nu = N21 = and A^u = /(z)A^22, where 



Recalling (4.5), we then find the following symmetries for the scattering data: 

A*(z*,r)=A*(l/z*,r), B\z*,T)^f{z)B*{\lz\T). (4.8) 

Note that N(z) is not invertible for z = ±x^^^,±X~^^'^- However, (4.8) is still valid at such 
values of z. Indeed, since <I>(z, solves (4.2), writing a Neumann series for €>(z, f) one finds 
^{+X^^^,t)\2 - 0, which implies that B*{±x~^^'^,t) =0. As a consequence, since A*(l/z*,r) 
and B*{\lz*,T) are analytic for z e C™, we can conclude that the limit as z — > ±X^^^ of the 
product /(z) B*(l /z*, 7") exists and is finite. 

The above properties now allow r*( 1 /z*) to be expressed in terms of the known functions, 
a{z) and b{z). For simpUcity, we consider the case in which no discrete spectrum is present. 
Consider first the case T - 00. The global relation in this case is simply given by (3.45). 
Replacing 1/z by z and using (4.8), we obtain 

a\z*)d{z) /(l/z)fc*(z'')^(z) „ 

A(l/z)= . B{llz) = -777- ze£>+in, 

A(l/z) A(l/z) 
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where 

A{z) = a{\lz)a*{\lz*)-vf{z)b{\lz)b*{\lz*). 
As a result, we can express the ratio R{z, T) =B*{\lz*, T) /A*(l Iz*, T) as 

R{z,T)^f{llz)^4ri zeD^in, (4.10) 

and we therefore obtain r*(l/z*) only in terms of known spectral functions. Now consider the 
case r < 00. The global relation in this case is (3.44). Replacing 1/z by z in (3.44) and using 
the symmetry (4.8) as before, we obtain 

a*iz*) a*iz*)A*iz*,T) 
We therefore see that the difference from the case r = oo is simply the appearance of an 
additional term in the RHS of (4.10). In Appendix E, however, we show that the second term 
in the RHS of (4.1 1) does not affect the solution of the IB VP for the IDNLS. Hence, even in 
the case T < oo, we can use (4.10) in the RHP (3.36). 



4.3. Discrete spectrum and soliton solutions 

Equations (3.35) imply that when the functions a(z) and d(z) possess zeros the matrices 
Mj7(z,f) are only meromorphic functions in D+ and D-, respectively. As a consequence, the 
RHP (3.36) formulated becomes singular. As in the IVP, however, it can be converted to a 
regular RHP by taking into account the appropriate residue relations. We assume that these 
discrete eigenvalues are all simple. More precisely, we assume that: 

• a{z) has simples zeros in Z)_in. We label such zeros ±Zj for 7 = 1,...,/; 

• d{z) has simple zeros in Z)+in. We label such zeros ±Aj for j =1, ... ,J'. 

We also assume that there are no zeros on the boundaries of these domains and that there are 
no common zeros of a(z) and d(z) in D+m . 

The fact that the zeros of a(z) and diz) always appear in opposite pairs is a trivial 
consequence of a(z) and d(z) both being even functions of z [cf. Appendix D]. Also, the 
symmetry Pn(t) = yqn(t) of the potentials implies that, corresponding to these zeros, there is 
an equal number of zeros of a*(l /z*) and d*(l /z*) in D+out and D_out, respectively, which we 
denote respectively by z, = 1/z* and A, - 1/A*. Thus, discrete eigenvalues in the IB VP can 
appear m two different kinds of quartets, namely, 

l+Zj, +-Zj]j^i , {+Aj, +Aj]jLi . 

Similarly to the IVP, from (3.27) and (3.36) we find the following residue relations; 

Res [M(-'^>] = a J M'^-'^'^Zj) , Res [Mi^'^>] = a, mI^'^^z,) , (4. 12fl) 

Res [M<^'^)] = dj M^^'^^KAj) , Res [M^-^^] = dj M^-'^KAj) , (4. 12^.) 

where 

aj = KjzY&-^''^^i^' , aj = Kjzf"e^''^^J^' , dj = KjA^Q-^''^^'^' , dj = KfAf^e^^^'^' , 

Kj = 1 mzj) b{zj)) , A; = vB*{Aj)l{a{Aj) d{Aj)) , Kj = i-z*)-^vKj , Aj = i-A*)-\Aj . 

and as customary Kj, Aj, Kj and Aj are referred to as norming constants. Note that since h(z) 
and B*(l/z*) are odd functions of z [cf. Appendix D], the norming constants Kj at z = ±Zj are 
identical, and the same follows for A^ at z = ±Aj. 



L 
+ 
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The RHP is now solved by removing the singularities, which is done by subtracting the 
residue contributions at the poles. As usual, the solution of the RHP then has additional terms 
compared to the case of no poles (3.37), and is given by 

M„(z,0 = I + / M:(f,0^^^d^,+ I (-!- Res[M-(z)] + ResLM^Cz)]) 

Im-'^ f-Z i=\\z-Zj z=Zj z-Zj z=Zj I 

- Z (-^ Res[M:(z)] + Res[M;(z)]], (4.13) 

where we defined Zj+j - -Zj for /' = 1,...,./ and Aj+j> - -Aj for /' = 1,...,/'. From the 
asymptotic expansion of (4.13) and the symmetries (D.8), we then obtain the reconstruction 
formula: 

qn(t) = -2 2 e-2''^(^^>'/:,M;^i ii(z;) - 2 2: e-2''^(^^>'A;M;^i ,^(Aj) + UO , (4.14) 

where q„{t) is given by (3.41). 

In the reflectionless case with v = - 1 , we obtain the soliton solution solving the following 
algebraic system of equations for M^^j ^^izj) and MjJ^^j u^-^jy- 

M-nto) = 1 + Z -aji-^ - -^)Kn(zj) + Z dji-^ ^Wn.iiC^j) 

j=i \zi-zi zi + Zjl j=i ^zi-Aj zi + Ajl 

Kxiizi) = i ajij^ + j^)K,nizj) + 2 dj(^ + ^)m;ii(A;), 

j=i \zi-Zj zi+Zjl j=i \zi-Aj zi + Ajl 

j=\ \Ai-zi Ai + Zjl j=\ \Ai-Aj Ai + Aji 

M-^a/)- i; a,(y^ + yl-)M-ii(z;)+ 2 djij^ + j^]Ml,,{Aj). 
j=i \Ai-Zj Ai + Zjl j=\ \Ai-Aj Ai + Ajl 

For a single quartet {+z\ , +Z\ }, the solution of the above system with 7=1 and J' =0 yields 
the one-soliton solution of the IDNLS as 

qnit) = e2«(«+i)^+2«;r+^] siiih(2Q;) sech[2((n +\)a-vt-S)], (4.15) 
where zi = e"'^'^ and 

w = cosh(2Q;) cos(2yS) - 1 , v = sinh(2Qr) sin(2j8) , 

1 . X 1 TT 1 

The soliton solution corresponding to a single quartet {±Ai, ±Ai} has an identical functional 
representation, which also coincides with the well-known one-soliton solution in the IVP. 

Note that the norming constants A/ contain the unknown scattering datum q-i(t) through 
the spectral functions A(z, f) and B(z,t). In general, this datum must be obtained by solving a 
nonlinear system of ODEs, as explained previously. In the case of linearizable BCs, however, 
Aj can be expressed only in terms of known scattering data. In particular, with T = oo, the 
global relation imphes 

A,- = f(\IAj)b\A*j)l[a{Aj)k{\IAj)\ , (4.16) 

where A(z) was defined in (4.9). This result can then be used in the residue relations (4.12Z?). 

Equation (4. 16) is a consequence of the fact that d{z} and A(l/z) have the same set of zeros in 
Z)+in, which in turn can be easily proved considering the analyticity of A*(l/z*) and B*(l/z*) 
with (4.9). 
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5. Continuum: linear and nonlinear Schrodinger equations 

In order to compare the solution of the IB VP in the discrete case to its continuum limit, and 
to appreciate the differences between the method for discrete problems and its continuum 
counterpart, here we briefly review the solution of IBVPs for the linear Schrodinger (LS) 
equation and the nonlinear Schrodinger (NLS) equation: 

.■^ + g-2v|<,|2^ = (5.1) 

(with v = 0,±l denoting respectively the Unear, defocusing and focusing cases), to which (1.1) 
and (1.2) reduce to in the Umit ^ -> 0. Note that, even though the IVP for (5.1) was solved 
in the early days of integrable systems for both vanishing [44] and nonzero [45] BCs, the 
IB VP on the half line was solved only recently [29]. Also, even though the IVP for the vector 
generalization of (5.1) was also solved early on in the case of vanishing BCs [34], the analogue 
problem with nonzero BCs was also only recently solved [38]. 

Linear Schrodinger equation: IVP and IBVP via Fourier methods. Consider first the initial 
value problem for the LS equation with x e R, ; > and q{x,Q) given. For simplicity we 
assume that q{x,Q) belongs to the Schwartz class, which we denote by 5(R). The IVP is 
trivially solved using the Fourier transform pair, defined as 

CO 1 oo 

q{k,t)= ^ Q-'''''q{x,t)dx, qix,t) = — J c'''mk,t)dk . (5.2) 

-oo -oo 

Use of (5.2) yields the solution of the IVP as 

1 oo 

qix,t) = — fe'^'"'-'^'^qik,0)dk. (5.3) 

Now consider the IBVP for the LS equation on the half line with Dirichlet BCs; i.e., x>0, 
t>0 and with qix,0) e 5(R"^) and qiO,t) € C(R"^) given. Employing the sine transform pair 

oo 2 ^ 

q^'\k,t)= fsm{kx)qix,t)dx, q(x,t) = - f sm{kx)q'~'\k,t)6k, (5.4) 

^ 

yields the solution of the IBVP as 

2 2 1 2 

q(x,t) = - f e-'^ 'smikx)q'''\k,0)dk+- [ e"'*^ ' sm(kx)g(k,t)dk, (5.5) 

where 



^0 ^0 



g(k,t) = 2ik /e'^'''^(0,f')d?' . (5.6) 



5.7. Linear Schrodinger equation: IVP and IBVPs via spectral methods 

An algorithmic method to obtain the Lax pair of linear PDEs was given in Ref. [25]. However, 
one can also obtain the Lax pair for the LS equation via the linear Umit of the Lax pair of 
the NLS equation, namely (A. 2). Let Q = 0(e) and take <S>(x,t,k) - \{x,t,k) to be a two- 
component vector. To leading order it is \{xj,k) - '^'^'Vo, where Vo = (yi,o,f2,o)' is 
an arbitrary constant vector. Choosing V2,o = 1 and substituting into the RHS of (A.2a) then 
yields the following equations for fx{x,t,k) = e'^'^^~^'^'hi{x,t,k) up to O(e^) terms: 

yU;f — ik'ii — q, + ik'^fi — iq^ — k'q , (5.7) 

where k' = 2k. One can now verify that enforcing the compatibihty of (5.7) yields the LS 
equation. Hereafter, for convenience, we will omit the primes. 
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Initial value problem. Introduce a modified eigenfunction il/{x,t,k) = e~'^''^~'^'^^{x,t,k), 
which satisfies the simplified Lax pair 

^, = , = e-'^'^'-'\iq, - kq) . 

It is then easy to obtain the solutions of (5.7) which decay as x ^ ±00 respectively as: 

X 00 

fx^^\x,t,k) = J e'*(*-^')^(/,Odx', tP-Kx,t,k) = -/ e'*^(^-^')^(x',f)dx'. 

-00 X 

(5.8) 

Note that i/^'^\x,t,k) are analytic for \mk ^ 0, respectively Also, on Im^ = it is 

H^^\x,t,k)-n^^\x,t,k) = e"'''q(k,t) = e'<'^^-*^''^§(/t,0), (5.9) 
where qik,t) is the Fourier transform of qixj): 

oo 

qik,t)= / e-*^'g(x',Odx'. (5.10) 

— CO 

Also, i^^'^\x,t,k) = 0{\lk) as oo in their respective half planes. Thus (5.9) defines a 
scalar RHP which is trivially solved via the standard Cauchy projectors over the real Une: 

Mx,f,fc)=2^ /e'(*'-^'^')^d;fe'. (5.11) 

Inserting (5.11) into (5.7) then yields (5.3) as the solution of the IVP. 

Initial-boundary value problems. We now consider the IB VP for the LS equation on the half 
line. Define simultaneous solutions of both the x-part and the f-part of the Lax pair: 

tP{x,t,k) = J e'*(^-^)-* <'-''\qix',t')dx' + iiq^ix',/)-kqix',t'))dt'] . 

(Xj.tj) 

In particular, consider the three eigenfunctions pf-'^{x,t,k), j = 1,2,3, defined by the choices 
(xi,?i) = (0,0), ix2,t2) = (00,0 and ix3,t3) = (0,7): 

H^^\x,t,k) = /e*(^-^^^(x',Odx' + /e*^-'^('-'')(j^;,(0,?') -A:^(0,O)d?' , 

(5.12a) 

OO 

tP\x,t,k) = -/ e*(^-^^g(x',Odx' , (5.12fc) 

X 

X T 

tP\x,t,k) = J e'^(^-^)^(x',Odx'-/e'^^-'^ '^'-''\iq,iO,t')-kqiO,t'))dt' . 

' (5.12c) 

Note that p'-^^ coincides with the eigenfunction in the IVP. As for p^^^ and p'-^\ they are 
entire functions of k. These eigenfunctions have the following domains of analyticity and 
boundedness: 

p'-^^: k€Cn, p'-^^: A;€Cm+iv, p''^^: keCi, (5.13) 
where Cm+iv is the lower-half plane. The two jumps on Im A: = and the jump on Re A: = 
(with Im A: > 0) then define a scalar RHP: 

p^^\x,t,k) -p^^\x,t,k) = e"'''-"'^' Fik,T) ReA = OAlmA>0, 

(5.14a) 

p'-^\x,t,k)-p'-^\x,t,k) = e'''''-''^'qik,0), lmk = A ReA<0, 

(5.14b) 

p'^^\x,t,k)-p^^\x,t,k) = e''''-'''^'iqik,0)-Fik,T)), lmk = A ReA>0, 

(5.14c) 
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1 OO -t 

q(x,t) = — / e*^-'*''^(ifc,0)dyt-— / e^''-'^'F(k,T)Ak. (5.17) 



where F(k, t) = ifi (k, t) - kfo(k, t), and with 

oo t 

q{k,t) = J e-%(x,Odx, fnikj) = Je'*'*' dlqix,t')\x=odt' . (5.15) 



The one-sided Fourier transform q(k, f) is analytic and bounded for Im < 0, while the 
transforms f„ik,t) of the boundary data are entire, and are bounded for Im k^ > 0. Moreover, 
qik, ?) ^ as A: — > 00 with Im A: < 0, and f„{z, ^ as A: — > oo with Im A^ < 0. The solution of 
the RHP defined by (5.14) is thus given by 

u.{x,t,k) = — I e""^ — -dA' I e'*^^ ^dA;'. 

IniJ^ k'-k 2ni^ k'-k 

(5.16) 

Inserting (5.16) into the first of (5.7) then yields the reconstruction formula: 

1 oo 1 

- / e*--'*''^(A,0)dA-- 

As in the discrete case, (5.17) still depends on the unknown boundary datum qx(0,t) via its 
transform in F{k,t). Integrating (5.7) from (0,0) to (0,7), (oo,r), (oo,0) and back yields the 
global relation as 

T CO CO 

Je''^'{iqMt)-kq(Q,t))dt + e''''^ J e-%(x,r)dx = / e-%(x,0)dx, 



which holds for Im A < A Im A^ < 0, i.e., A e Cm- In terms of the spectral data: 

ifi (A, T) - kfoik, T) + e'*^''§(A, T) = qik, 0) , VA € Cm ■ (5.1 8a) 

Using the the transformation A ^ -A, which leaves fn{k,t) invariant, from (5.18a) we obtain 
ifiik,T)+kfoik,T) + e''^^qi-k,T) = qi-k,0) VAeQ. (5.18fo) 

We then solve for f\(k, T) and insert the result in (5.17). [The first term in the RHS of (5.18fc) 
yields a zero contribution to the solution.] Thus, the solution of the IB VP is given by 

1 oo 1 

q(x,t) = — re'^^-'^''$(A,0)dA-— [ e'*^^-''=''[§(-A,0)-2A/o(A,r)]dA. 
2^ -co 2;r^j 

(5.19) 

Note that one can replace /o(A, T) with /o(A, t). Also, the second integrand in (5.19) is analytic 
and bounded for Im A > A Im A^ < 0. Thus, one can deform the integration contour on 
the second integral onto the real A-axis and recover the sine transform solution (5.5). Unlike 
sine/cosine transform approaches, however, the present method can be applied to solve IB VPs 
with more compUcated BCs, as we show next. 

Robin BCs. Consider the IB VP for LS equation with Robin BCs: 

^^(O, t) + qAO, t) = hit) , (5 .20) 

with hit) given and where a e C is a nonzero but otherwise arbitrary constant. In a similar 
way as shown in Appendix C for the discrete case, one obtains [23, 25] 

Fik,t) = plIl + P±,^^''qi-k,t), (5.21a) 
k — ia A - ia 

where 

Gik,t) = 2ikhik,t) - ik+ia)qi-k,0) (5.21fo) 
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contains the known portion of (5.21a) and where hikyt) is defined according to (5.15). Then, 
again following similar steps as in the discrete case, one obtains the solution of the IB VP as: 

q{x, t)=^l e'*^-'*''§(A:, 0) dA: - f e*^"'*"' dk + ivae-""'^'"''Giia, t) , (5 .22) 

2n _^ Ln k-ia 

where Vq, = 1 for -7r/2 < argor < 0, v^^ = 1/2 for arga = 0, -7r/2 and Vq, = for < arga < 37r/2, 
and where the integral along d£.\ is to be taken in the principal value sense when arga = 
0, -7r/2. (The last term in the RHS of (5.22) is missing in Refs. [23, 25]. One can easily show, 
however, that without this term q(x,t) does not satisfy the BC at x = 0.) 

5.2. Nonlinear Schrddinger equation: initial value problem 

As in the linear case we assume that q(x,Q) € S(R.). Recall that the Lax pair for the NLS 
equation (5.1) is given by (A.2) with p{x,t) = vq*{x,t). For the present purposes, we consider 
^(x, f, ^) to be a 2 X 2 matrix. 

Analyticity. Introduce a modified eigenfunction which has a well-defined limit as x —> +oo: 
^i(x,t,k) = 0(x,f,yt)e-'"("'''*''"3 , (5.23) 

with 8(x, t, k) - kx- 2k^t. Note /i(x, f , k) satisfies the following modified Lax pair: 

Hx-ik[a-3,ii\=Qn, jut + 2ik^ [0-3 , //] = H// . (5 .24) 

Then, letting iu{x,t,k) = e'^*3T(x,r,A:), we obtain the simpUfied Lax pair: = e"'^^3(Q)*P 
and *l'r = e"'*''3(H)*i'. We then define the lost eigenfunctions as the solutions of (5.24) that 
reduce to the identity as x — > ±00: 

X 

tf^\x,t,k) = I + / e''^''-'^^^\Q{x' ,t)n^^\x' ,t,k))Ax' , (5.25a) 

—00 

CO 

fi^^\x,t,k) = I -/ e'''^''-'">^^iQ(x',t)i^^^\x\t,k))dx'. i5.25b) 

X 

We have the following regions of analyticity and boundedness [6]: 

where i/^\x,t,k) - (/z^^'^' as before. The analyticity properties of <i>^i\x,t,k) = 
H^'>ix,t,k)e'^"'\ 7 = 1,2, follow trivially. 

Scattering matrix. Note detO^^^ = det/^^^^ = 1 for j - 1,2. Thus O^'^ and O^^^ are both 
fundamental solutions of (A.2) Vk e R. Hence 0^^\x,t,k) = 0^'^\x,t,k)A{k), where A(k) is 
the scattering matrix. Equivalently, 

^i''^\x,t,k) = ii^\x,t,k)e^''^kik). (5.26) 

Note that A(^) is indeed independent of time, and det A(A:) = 1 . Moreover, 

oo 

A(A;) = I + / e-'^''''-^'^'^^\Q{x,t)n^^\x,t,k))Ax, (5.27) 

— oo 

and 

aii(A:) = Wr(<l)*'"^\(D^2'^^), fli2(/t) = Wr((D^''^^a)<2'*^), (5.28fl) 

fl2i(;t) = -Wr(0<'"^>,(I)<2'^'), fl22('t) = -Wr((D(''^\0<2"^'). (5.28/7) 

Thus, an{k) and a22{k) can be analytically continued respectively on ImA: < and ImA: > 0, 
but ai2{k) and aiiik) are nowhere analytic, in general. 
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Symmetries. When p{x,t) = vq*{x,t), with v = ±1, the scattering problem (A.2) admits an 
involution expressed via the matrix cTy in (A.lla): if 0{x,t,k) is a solution of (A.2a), so is 

<D'(x,f,yt) = o-^(D*(x,f,F). (5.29) 

Comparing the behavior of the Jost eigenfunctions as x — > +oo we then have 

Q)^j-^\x,t,k) = (Tyi^^j-'^\x,t,k*)T , ^^^•'^\x,t,k) = V(Tyi^^j-^\x,t,k*)T , 

(5.30) 

for j =1,2. Hence the following relations hold for the elements of the scattering matrix A{ky. 

a22ik) = al^ik*), a2iik) = val2ik*). (5.31) 
Note that, since detA(A:) = 1, (5.31) imply |aii(A:)p - v|ai2(A:)p = 1 VA; e R. 

Asymptotics. The asymptotics of the Jost solutions as A ^ oo in their half planes is: 

1 1 ^ 

/i<'\x,f,)fc) = I-— 0-3Q+ — 0-3 r q(x',t)p(x',t)dx' + 0(l/k^), (5.32a) 
2ik 2ik 

1 1 00 

lP\x,t,k) = I - — 0-3Q- ttc^3 r q{x',t)p{x',t)dx' + Oil/k^). (5.32fc) 

2ik 2ik •'^ 

Moreover, from (5.28) and (5.32) one also obtains 

CO 

a22ik)=l r f qix,t)pix,t)dx + Oil/k^). (5.33) 

'^ik -•'00 

Inverse problem. The inverse problem is the RHP defined by (5.26) for A: e R: 

U-{x,t,k) = M+(x,/,yfc)(l - J(yfc,/)), (5.34) 
where the matrix- valued sectionally meromorphic functions are 

M-(x,a) = (/.(^•^>(x,a),^^^^^^%^), M-(x,a) = (^^^^^^%^,ju(2'«'(x,a)), 

V a22{k) I \ an{k) I 

the jump matrix is 

^^^''^"i-e-2'%(fc) j' 
and the reflection coefficients, defined VA e R, are 

p\{k) = a2\{k)lan{k), p2(k) = anik) I a22{k) . 

Of course (5.31) imply pi{k) = vp*2{k*) when p{x,t) = vq*(x,t). In the absence of a discrete 
spectrum [i.e., if aii(A) ^ VlmA; < and fl22(A) VlmA > 0] the matrix functions 
M-(x, t,k)-\ are sectionally analytic in their respective half planes, and they vanish as A: ^ 00. 
Therefore the RHP (5.34) is solved via the Cauchy projectors P^, as for the linear case: 

1 JC/k' t^ 

U^{x,t,k) = \ + — r U^{x,t,k')-i-^dk' . (5.35) 
2m 4x, k' -k 

The asymptotic behavior of M(x,;,A) as A —> 00 is easily obtained from (5.35): for ImA: > 0, 

1 °° 

M^{x,t,k) = \ r M+(x,?,A:')J(A',?)dA:'-i-0(l/A^). (5.36) 

Comparing the (l,2)-components of (5.36) and (5.32) then yields the reconstruction formula: 

1 CO 

q{x,t) = - / e2'(^^-2^''V2(A)(/i<^\x,?,A))iidA. (5.37) 



IBVPsfor discrete evolution equations: DLS and IDNLS 29 
Linear limit. If Q{x,t) = 0{e) one has n{x,t,k) = I + 0(e) and, to 0(e), 

oo 

A(A:) = I+ / e-'('=^-2'=''''^3Q(x,0dx. 

— CO 

From here and (5.37) one then obtains, to 0(6:), 

■J oo CO 

q{x,t)=- r e^'^''''-^'^'^p2{k)Ak, p2{k)= r e-2%(x,0)dx, 

" — oo — oo 

which, with the familiar rescaling k' - 2k, coincide with the Fourier transform pair (5.2). 
5.3. Nonlinear Schrddinger equation: initial-boundary value problem 

We now discuss the IB VP for the NLS equation (5.1) on the half Une. As in the linear case, 
we assume q(x,0) € .S(R+) and ^(0,0 g C(R+). 

Eigenfunctions and analyticity. Introduce three Jost eigenfunctions as the solutions of (5 .24) 
that reduce to the identity respectively at {x,t) = (0,0), {x,t) — > (ooj) and (xj) = (0,7): 



tf\x,t,k) = \+]e''^''-''^^^{Q{xf,t)i/^\x',t,k))Ax' 

.jgi[*A;-2fe2a-O]^3(H(0,?',A:)//(i'(0,?',fe))dr', (5.38a) 





+ Je" 





tP\x,t,k) = \-Je'''^''-'^^^^(Qix',t)fx'-^\x',t,k))dx', (5.38fo) 

X 

H^^\x,t,k) = \+] e''^'-'^'^^\Q{x' ,t)n^^\x' ,t,k))dx' 

T 

_ jre«*^-2^('-'')]'t3(H(0/,A:)/i(3)(0,;',A:))dr' . (5.38c) 

t 

Note that n^'^\x,t,k) and i/^\x,t,k) are entire functions of k, while iP'\x,t,k) coincides 
with (5.25fo). Moreover, (5.38) imply the following domains of analyticity and boundedness: 

/i(i'^>: Cm, //^i-"': Cn, tP'"^^ : dv, tP'"^ : d, 

Scattering matrices. We still have det<l>*^\x,f,fc) - 1 for all x,t€ R"*" and for all j - 1,2,3. 
Hence the matrices O*-'' (x, t,k), j - 1 , 2, 3 are three fundamental solutions of the Lax pair (A.2), 
and they must be proportional to each other. In terms of the modified eigenfunctions: 

//(2\x,f,yt) =//(^'(x,f,fc)e'(*"-2*'"*3s(fe), (5.39fl) 

lJ^^\x,t,k) = //("(x,f,yt)e'('^^-2*'"*3S(/t,r). (5.39^) 

Note that the first column of (5.39a) is defined VA: € Q+n, the second column Sk g Cm+iv 
and (5.39fo) holds VA; g C. Also, dets(A:) = detS(fe,r) = 1. The scattering matrices s{k) 
and S{k, T) are obtained from the boundary values of the eigenfunctions, namely, VA: g C, 

s{k)=tP\Q,Q,k), S(A:,T) = (e2'*^^%(i>(0,T,fe))"V (5.40) 
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Then, from (5.38) we have the following integral representations of the scattering matrices: 

sik) = I - / e-'*^'^3(Q(x,0)/|(2^(x,0,yt))dx, (5.41a) 



T 

S-\k,T) = I + / e^'* '^^(H(0,t,k)n^^\0,t,k))dt. (5Alb) 



These imply that: Sl(^) and SR(k) are analytic respectively for k e Ci+n and k e Cui+iv, and 
their restriction to these domains are continuous and bounded on the boundary; S{k,T) is 
entire, and S^ik, T) and Sfi(A:, T) are bounded respectively for k e Cn+iy and k e Ci+m- 



Symmetries, discrete spectrum and asymptotics. When pix,t) = vq*ix,t), (5.29) still holds, 
as does (5 .30) for 7 = 1,2,3. This implies that the scattering matrices can be expressed as 



sik) = 



a(k) vb*(k' 



1\ 

■) r 



Sik,T) = 



A{k,T) vB*(k 



\T)\ 



b{k) a*ik*) I' ' ' ' \Bik,T) A*{k* 

The properties of a{k), b{k), A{k, T) and Bik, T) follow trivially from those of s(A:) and S(A:, T). 
Also, one can show that f/j\x,t,k) = I + 0(l/k) for j = 1,2,3 as — > oo in the respective 
domains of boundedness of their columns. The asymptotics of the eigenfunctions then 
determines that of the scattering matrices. In particular, aik) = 1 + 0{l/k) and b{k) = 0(1 /k) 
as A: ^ 00 in Ci+n, and Aik, T) = 1 + 0(1 /k) and Bik, T) = Oil/k)ask-> coin Cm+iv- 



Riemann-Hilbert problem, solution and reconstruction formula. Equations (5.39) allow us 
to formulate the following RHP: 

M-ix,t,k) = M^ix,t,k)i\-Jik,t)), k€L, (5.42) 

with L = 5Ci U dCui ^ L1UL2UL3UL4, where 

Li=CinCn, L2 = CnnCni, L3 = CmnCiv, L4 = CinCiv, 

and where 



M^ix,t,k) ■ 



,(2,Z.) 



0,R) 



' dik) / ' 



A:€Ci, 



M ix,t,k) ■■ 



[ \a*ik*y' 

The jump matrices Jjik,t), each defined for k e Lj, are: 



,(2,L) A* 



,(3,L) 



' aik) 



,(2,R) 



d*ik*y'' 



fceCn, 
A; € Civ . 



hik,t) ■ 



ve2'»r*(A*) 




-e-2'«r()fc) 



^2ik,t) -y_^^.2U)^,^^^ ^|^(^)|2 j. 

J4(A:,0=I-(I-Ji)(l-J2)"kl-J3), 



and the reflection coefficients are 

vb*ik) 



yik). 



dik) = aik)A*ik*, T) - vbik)B*ik*, T) , Yik) ■■ 



Bik,T) 



aik) a*ik*)d*ik*) 

Note that dik) is defined VA e Cj+u, Tik) for A e L3 U L4 and yik) Vfc E R. Their asymptotics 
as A: ^ 00 follow trivially from those of sik) and S(A:, T). As a result, M(x,f,A;) — > I as A — > 00. 
Hence, in the absence of a discrete spectrum [that is, assuming that aik) and dik) have no zero 
respectively for A e Cn and A e Q], the RHP (5.42) is solved by Cauchy projectors: 



1 t\ 
U^ix,t,k) = \ + — \U^ix,t,k')-i-^6k' . 

Ini ■£ A:' - A: 



(5.43) 
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Substituting the asymptotic expansion for M{x,t,k) into the x-part of the Lax pair and 
comparing the (1,2) components, we have 

q{x,t) = -2i Urn k(M(x,t,k) - I)j2 . (5.44) 

Using the asymptotic expansion for W\(x,t,k) as — > oo, from (5.43) and comparing the (1,2) 
components, we obtain the solution of the IB VP for the NLS equation as 



qix,t) = ^ / ye2'^(^''''^')r*(yt'*)M|i(x,f,r)dA:' 

— / e^'^^'''''*'V(A:')Mti(JC,?,A:')dA;'--/v|r(-A;')pM+2(jc,r,-A:')dA:'. (5.45) 



Linear limit. Supppose that Q - 0(e). From (5.38) and (5.45) we have ji -\ + 0(e) and 
M = I + 0(s). Also, (5.41) imply y(k) = -q(2k,Q) + O(e^) and d(k) = 1 + 0(e\ as well as 

r*(^*) = v(2^/o(2^, T) - if (2k, T)) + 0(e^) . 
Thus (5.45) yields, to 0(e), 

1 1 OO 

q(x,t)^ - [ e2'(*'^-2*"'>(2/fc7o(2/t',r)-//i(2fc',r))dfc' + - re2'(*'^-2*"^)^(2/t',0)d)t'. 

Performing the change of variable 2k' = k, we then see that, to leading order, this expression 
yields exactly the solution of the linear Schrodinger equation on the half Une, namely (5.17). 

Global relation and Dirichlet-to-Neumann map. Equations (5.41) involve all initial and 
boundary data for Q(x,t). These values are not all independent, however, since they satisfy 
the global relation 

T CO 

/ e^^'^-^s (H(0, t, k)n(0, t, k)) dt + e^'^'^^s J ^-ikxe-^ ^q^^^ j^^^^^^ j^y^ 

CO 

= / er'^"^^ (Q(x, 0)n(x, 0, k)) dx . (5 .46) 



When (5.46) is evaluated with ji = j/^\x,t,k), its first column is defined VA: € Ci+u, its second 
column VA: e Cm+iv. Moreover, when n(x,t,k) - i/^\x,t,k), the RHS of (5.46) equals I -s(A:). 
Using (5.39fc) in the LHS, one then obtains a relation between the scattering matrices: 

S-^(k, T)s(k) = I - s^''^^^^G(k, T) , (5.47) 

where 

oo 

G(k,t) = J e-''"'^^{Q(x,t)tP\x,t,k))dx, 



and Gi,(k,t) and G^(fc,f) are analytic respectively for k e Ci+u and k e Cm+iv, and continuous 
and bounded on the boundary of these domains. In particular, for k e Cm+iv we have 

CO 

A*(k* J)b*(k*)-B*(k* J)a*(k) = -ve'^'''^'^ fe-^'''''q(x,T)nf^(x,T,k)dx. 



Since the integral term in the RHS is of 0(1 /k), as A: ^ oo in Cm, integrating along dCm we 
obtain the following integral relation: 

/ ks-'^'''^'{B*(k*)-r(k*)A*(k*))dk = 0, (5.48) 

dCiu 

where r(k) - b(k)/a(k). As shown in [28], this relation can be solved to obtain the Dirichlet-to- 
Neumann map, which expresses the unknown boundary datum qx(0,t) in terms of the known 
one, q(0,t). 
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Linearizable BCs and soliton solutions. One can write S{k,t) = 3) ^{k,t)e ^''^"^3"'3^ where 
^ik,t) = O''^\0,t,k) solves the f-part of the Lax pair {A.2b) for jc = 0, namely 

^t + 2ik^cr3^ = H(0,t,k)^, (5.49) 

with 0(k, 0) = I . The matrix ^(-k, f) solves an equation identical to (5 .49) except that H (0, t, k) 
is replaced by H(0,f, —k). If there is an invertible time-independent matrix H{k) such that 

N(A;) (2jA^o-3 - H(0, t, k)) = {lik^cr^ - H(0, /, -k)) H{k) , (5.50) 

it then is easy to see that ^{-k, t) = N{k) 3>(A:, t) (k) . One can show that a suitable matrix 
N{k) only exists for homogeneous Robin BCs, namely, 

with ;t' £ arbitrary. In that case, (5.50) implies N\2 - Nix = and A'n = f(k)N22, where 
f(k) = -(2ik-x)/(2ik+x), which in turn imply A*(k*,T) = A*(-k*,T) and B*(k*,T) = 
f(k)B*(-k*, T). From here, similar arguments to those used in the discrete problem can be 
apphed to the analysis of linearizable BCs. 

As in the discrete case, the poles for the IBVP occur at the zeros of a{k) in Cu and those 
of d{k) in Ci, plus their complex conjugates in Cm and Civ [29]. Each of these pairs of zeros, 
by itself, generates the well-known one-soliton solution of NLS: 

q{x, t) = 277e2'f^-4'(f'-"')'+'('^-'^/2) sech(277x - 8f //f - 26) , (5.51) 

where ki-^ + irj is the zero of a{k) or of d{k), and C\ = Irje^^'^"^ is the norming constant (see 
[29] for further details). 

6. Conclusion 

In conclusion, we have demonstrated a method to solve initial-boundary value problems for 
linear and integrable nonlinear discrete evolution equations. We have done so by solving the 
IBVP for the discrete linear Schrodinger (DLS) and integrable discrete nonlinear Schrodinger 
(IDNLS) equations on the natural numbers. Moreover, we have illustrated the similarities and 
differences between the method for differential-difference equations and PDEs by showing 
explicitly the correspondence between the discrete and its continuum limit. While the 
differential form representation of the continuum is lost, the essential ideas of the method 
can be carried over to the discrete, but the actual implementation of the method presents some 
additional difficulties. In particular, the jump location in the nonlinear case differs because of 
the rescaling z' - in the dispersion relation a)(z) when going from the linear to the nonlinear 
case. This is a significant difference from continuum limit, where the jumps in the nonhnear 
case are given by the union of those for the linear problem and its adjoint (cf. sections 3 
and 5). Also, the limit —> oo in the continumm becomes z — » (for ImA: > 0) and z — > oo (for 
Imk < 0) in the discrete. As a consequence, the behavior of the eigenfuncions and spectral 
data as z — > in the discrete problem must also be studied in addition to that as z — > oo. This is 
the why the point z = plays such a special role in the discrete problem, similarly to Ref. [1], 
and is one of the reason why discrete problems are more complicated than their continuum 
counterparts. 

For the DLS, in addition to solving the IBVP with Dirichlet-type BCs we have shown 
that, contrary to Fourier series approaches, the method can deal with more complicated kinds 
of BCs just as effectively. For the IDNLS, in addition to solving the IBVP (showing explicitly 
how to eliminate the unknown boundary datum), we have characterized the linear limit, the 
linearizable BCs (showing how they fit within the 1ST framework), and we have obtained 
explicitly the soliton solutions. 
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It should be clear that, similarly to the continuum, the method can be generalized to solve 
IB VPs for both the DLS and IDNLS equation defined on a finite set of integers. It would also 
be straightforward to generalize this method to any discrete hnear evolution equation and to 
other integrable discrete nonlinear evolution equations. 

Several interesting questions can now be effectively addressed using the present method. 
For example, one can use the expression for the solution to study its long-time asymptotics, 
using the Deift-Zhou method [17], or to study the "small dispersion" or "anti-continuum" 
limit (i.e., the limit h — > oo), e.g., using the Deift-Venakides-Zhou method [16]. Doing so is a 
nontrivial task, however, which is beyond the scope of this work. 
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Appendix A. Notation and frequently used formulae 

We denote the closure, interior and boundary of a domain D respectively by D, D" and 
dD, where as usual dD is oriented so as to leave D to its left. We also occasionally refer 
to punctured regions of the complex plane, which we denote as R^'^^ — R-{Q]. As usual, 
[A, B] = AB - BA is the commutator of two matrices A and B. We use a superscript asterisk 
to denote the complex conjugate z* of a complex number z, and |zp = z*z. Throughout, 
R+ = {xeR: x>0) andRj =R+U{0). Similarly, N = {1,2,3,...} andNo =NU{0). Finally, 
we denote by Q, . . . , Qv the first, second, third and fourth quadrants of the complex plane: 
Ci = (A: e C : ReA: > A Imk > 0}, etc. Similarly, we denote by Q+n = (A e C : \mk > 0} and 
Cm+iv — {k eC'.lmk <Q] the upper-half and lower-half planes, respectively. 

The nonlinear Schrodinger (NLS) equation (5.1) is a reduction of the system 

iqt + qxx + ^q^p ^0, (A. la) 

-ipt + Pxx + 2p^q = 0. (A.lb) 

That is, (5.1) follows by imposing pix,0) = vq*{x,0) in (A.l), which then implies that 
pix,t) = vq*{x,t) > and qix,t) is a solution of (5.1). A Lax pair for (A.l) is given by: 

a);c-(lo-3<D = Qa), (Ala) 
Of -I- 2ik^cr30 = H <D , (A.lb) 
where ^{x, t,k) is either a 2-component vector or a 2 x 2 matrix, and where 

-3=(i Oix,t)^[l l], (A.3a) 

Hix,t,k) = .-.3(0.-0^2.0 = ( (^-3^) 

(The present pair differs from that in Ref. [26] by the rescaling k —> -k, and from that in 
Ref. [6] by A: — > -A and t — > -f.) Similarly, the integrable discrete NLS equation (1.2) is a 
reduction of the system of differential-difference equations 

iqn + iqn+l - 2qn + ^n-l) - qnPn{qn+l "H ^„-i) = , (A.4a) 
iPn + {Pn+l-'^Pn + Pn-l)- Pnqn{Pn+l + Pn-\) = 0. {AAb) 
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That is, imposing p„(0) = vg*(0) on (A.4) yields p„{t) = vql(t) "it > 0, with q„{t) 
satisfying (1.2). In the literature, the name Ablowitz-Ladik (AL) is associated to both (1.2) and 
(A.4). To avoid confusion, here we will simply refer to (1.2) as the IDNLS equation, reserving 
the name AL for the more general system (A.4). A Lax pair for the AL system (A.4) is: 

a>„+i-za)„ = Q„a)„, (A.Sa) 

a>« - i(z - l/z)V3 <Dn = H„ <D„ , (A.5fc) 
where <I>n(z, t) is either a two-component column vector or a 2 x 2 matrix, and where 

(AM) 

In sections 3 and 5 we make frequent use of the integrating factors 

Z(A) = ZAZ-i=f "'"^'l ^3A = ("" -""'A, (A.7) 

\fl2l/z fl22 / \ --221 022 / 

e''"^3(A) = e'e'-3Ae-'«'-3=( «ii e^'^a^ \ 
\ e ^'''aai 022 / 

For any matrix A, we write A = (A^^\ A'^^^), where the superscripts L and (Icfl and right) 
denote respectively the first and second column of A. We also write A = A^ + Ao, where A^ 
and Ao denote respectively the diagonal and ofl'-diagonal part of A. Note that 

(Aiu)d = Aoyuo + Aoi^o , (Aju)o = AoiUD + Aduo , (A.9a) 

(Q/z)d = Quo , (Q/i)o = QHD , iA.9b) 

and in particular 

H„^Diz,t) = -io-3Q„Q„-i, H„,o(2,0 = j(Zo-3Q„ + Q„_iZo-3). (A.IO) 

Note also that ZAo = A^Z"' and cr3Ao - -Aqctj. 

The "involution symmetry" of the scattering problems of NLS and IDNLS is expressed 
through the matrix 

cry = l^_ I]. (A.lla) 



^ V 0^ 

That is, when pix,t) = vq*(x,t) in (A.3a), or p„(t) = vq'^it) in (A.6a), it is, respectively: 

CTyQ* = QcTy , CTvQl = QnCTv • (A.llb) 

Note also that CyZ = Z~^ay, CTyaj, = -cr-jo-y, and cr~^ =cr'y = very. 

When discussing the asymptotic behavior of the eigenfunctions, the behavior of the 

matrix product AZ motivates the following definitions: for any matrix A = (A*^^', A*^^^), we 
write A = 0(1'") as z (0, 00) if A<^) = 0(2™) as z ^ and A^^) = 0(1 /z'") as z ^ 00. Similarly, 
we write A = 0(Z'") as z ^ (00, 0) if A^^^ = 0(z"') as z ^ ex? and A^^^ = 0(1 /z") as z ^ 0. 

Appendix B. Spectral analysis of the r-part of the Lax pair of the DLS 



The inversion formulae for the spectral functions (2.14) in the linear problem can be obtained 
by performing spectral anlaysis of the individual parts of the Lax pair (2.6). The first of (2.17) 
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can be derived from similar steps as in section 2.2. As for the second of (2.17), consider the 
following spectral problem 

lit + ioi{z)n = m, (B.l) 

where cj{z) = 2- (z+l/z). The Jost solutions are easily obtained, and are: 

fi^'\z,t) = /e-'"(^'('-^>/(/)dr', ^^^\z,t) = -/e-'"(^>(^-''>/(?')d?'. 

t 

Note that /x^^^ and ju^'^^ are analytic for ziD+ and z e D+, respectively, where D+ is the same 
as in section 2.3. Also, the jump condition is 

ju^i^ - n^^^ = Q-^'^'fiz, T), z€dD+, (B .2) 

where 

fiz,t) = j^''^'^''fit')dt'. 



Using integration by parts, one can show that fi^ - Oil/z) as z oo in their corresponding 
domains. Hence the solution of the RHP (B.2) is given by 

Substituting this into (B.l), we then find the reconstruction formula 

Recall that dD+ - dD+i^ U dD+out- Also note that dD+ia can be deformed to <9£>+out by letting 
z — » 1/z, and (o{z) and fiz,t) are invariant under this transformation. After some algebra, we 
then obtain 

/« = ^ / (4-i)e""^'VUr)dz. 

Replacing fiz,t) by qnit), we finally obtain the second of (2.17). 

Both of (2.17) could also be obtained by more direct methods. The first of (2.17) of 
course just defines the coefficients of the principal part in the Laurent expansion of qiz,t). 
As for the second of (2.17), it can be obtained as follows. Define q(t) to be the function 
which equals q„(t) for <t <T and is otherwise. Also, let 2(w) = J _^e"^'q„(t)dt be its 

Fourier transform. Then, for all < f < J it is qn(t) = (l/2n) J'^^&~"^*Q{aj)(kj. Note however 
that the transformation z oj(z) maps dD+out onto the real a>-axis, with (o{z) decreasing 
monotonically as Rez increases. Moreover, Q(a>(z)) = f„(z,T). Hence we can rewrite the 
previous integral as qn(t) = (l/2n) oj'(z)&~''^^'^'Mz, T)dz. 

•J OU+nnt 



Appendix C. IB VPs for DLS with Robin-type boundary conditions 

Consider the DLS equation (1.1) forneNo and/eR"^ with mixed BCs. The spectral transform 
of (2.27) yields, VzeCt^l, 

f-iiz,t)-aMz,t) = kz,t), (C.l) 

where the fj(z,t) are given by (2.14), and h(z,t) is defined similarly. Recall that the 
reconstruction formula (2. 19) contains the quantity F(z, ?) = iizfoizj) - f-iiz,t)). Use of (C.l) 
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and the transformed global relation (2.23) allows one to eliminate fo{z,t) and f-i{z,t) and 
express F{z,t), for all < |z| < 1, as 

G(z,t) z- 



F(z,t) ■ 



q{llz,t). 



(C.2) 



Ijz-a l/z-a 

where Giz,t), which contains the known portion of the RHS, was given in (2.29). Now 
recall that, in (2.19), Fiz,t) is integrated along dD+in. Three possible situations can arise: 
(i) a e D+out, (ii) a e f^-D+out, (iii) a i 5+out- We discuss each of these cases in turn. 

If a ^ D+ouu the denominator of (C.2) never vanishes in D+in. Thus the second part of the 
RHS of (C.2), when inserted in (2.19), gives rise to an integrand that is analytic and bounded 
in D+in. Hence, that part of the integral is zero. As a result, the solution of the IB VP is simply 



qnit)-:^. § z«-ie--(^)'§(z,0)dz--L / z«-ie--(* ^J^dz, 

\Z.\—l f'i-'+in 



(C.3) 



with G(z,t) again given by (2.29). Now suppose a g D+out- In this case l/z-a vanishes at 
z = l/or e D+in. Even though each of the two terms in the RHS of (C.2) has a simple pole at 
this point, however, their sum is finite there, since PizJ) is analytic in C^^^ Thus, 



55 f 



z-a 
l/z-a 



q{l/z,t)dz= Res 



z=l/a 



z"-'^^q(l/z,t) 



l/z-a" 



= Res 

z=l/a 



£1-1 ^-i(u(z)t 



G(z,t) 



l/z-a 



-n-l -i(i)(a)t 



G(l/a,t), 



which imphes the solution of the IB VP as 

^„(;) = -L (f 2«-ie-'<^»§(z,0)dz--^- f z"-ie-'"(*^^dz-a-"-ie-'"('^)'G(l/a,0. 
2m, r_, 2m^;i l/z-a 

(C.4) 



Finally, if a € 5D+out. the pole is along the integration contour. In this case one should go 
back to the RHP and subtract the pole contribution. In this way, the solution of the IBVP can 
be obtained as 

qnit) = 

2m 



£-1 Q-i(o(z)t ^ 



lzl=l 



q{z,mz - r z"-ie-''"«' dz - ia-"-ie-'"(«)'G(l/a,0 . 

2m„;L 1 Z-a 2 

(C.5) 



Combining (C.3), (C.4) and (C.5) one then obtains (2.28). 



Appendix D. Asymptotic behavior of the eigenfunctions of the IBVP 

DLS. We first compute the asymptotics for for n = (where no summation is present), then 
consider the casen > 1. Note that a>(z) = -l/z + 0(l/z^) as z ^ 0. Integration by parts yields, 
as z — > with Imz < 0, 

cf,^^\z, t)^q-i (t) - e-''^'^'q-m + Oiz) , 
while as z — > with Imz > it is 

cf>^^\z,t) = <7_i(0-e''^«('-^)^_i(r) + O(z). 

Using these in (2.13) with n > 1 we then have immediately (p^J\z,t) = qn-iit) + Oiz) as 
z ^ with Imz < for 7 = 1 and Imz > for j = 3. Note also that 0^ \z,?) - 0[f\z,?) = 
_g-M*(^_j(0)_e'"(^>^^_i(r)) + C>(z) as z ^ 0, implying that the ratio F(z,T)/z in (2.14/7) 

(2) 

remains bounded as z ^ along the real axis. As for <p„ iz,t), (2.9) implies immediately 
(f>^^\z,t) = Oil/z)a&z^oo. 
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IDNLS. The determination of the asymptotic behavior in the nonlinear case is considerably 
more involved, and requires the use of a Neumann series approach: 

f4P(z,t)= ifi';!'^\z,t). (D.i) 

m=0 

We now show that, Vn e No, m>0 and j = 1 , 3, as z ^ (oo, 0) it is 

fi%'"-'\z,t) - 0(Z-2'«), //g""''(z,0 = OCZ-^™-'), (D.2a) 
Hf^\z,t) = 0(2-2-), ^,^^\z,t) = 0(Z-2™-i). (D.2Z.) 

The proof proceeds by induction. Consider ui'\z,f) first. Separating (3.25a) into its diagonal 
and off-diagonal components then yields jl^ iz,t) = I and fr^'gizj) = O, as well as 

n-l 



/il'r''fcO= 2 Q„KOA'lVr'^feOZ-'+/(Ho,D/i[,^o"'' + Ho,oA*[,'o"'"^^)(z,?')d^', (D.3a) 

n'=0 ' ' ' 

^'^•;''\z,t) - Q„'(f)//l^S^(z,0Z-2("-"''-i 

n'=0 

j^-Mz)(«')^3(Ho,o/i<y + Ho,z,/i[,y )(2,^')Z-2«d?' . (D.3^) 







Note that 



^ I = -Z-2 + 0(Z-'*), asz^(oo,0). 



2w(z) 

First consider the case n = 0. Using integration by parts in (D.3^), we obtain, as z — > ('>o,0), 
4'o"'\z,t) = {Q-^(t)^^'^■;;\z,t) - ^-''^'^'^^ [Q_i(0);.[,;^'"(z,0)]}Z-i 

+ {{QoQ-i)m'^^'^\z,t) - e--(^)'-3[(Q,Q_j)(0);,(i^"')(^,0)]}Z-2 , (D.4a) 
plus higher order terms. Substituting (D.4a) into (D.3a) with n = 0, one finds 

/'oD^^fcO = -j7cr3(QoQ-i)(?')/^[,'a"^(z,?')d?' + i7cr3Qo(f')/^[,'r'^fc^')Zdr' 



- j7cr3Q-i(?');«[,'r'^fc^')Z"' d?' . (D.4fo) 



Using (D.4) one can then obtain (D.2) for n-0 and all m e No inductively. Note also that, for 
m = 0, (D.4a) yields (3.32^). Similarly, repeating the same arguments, one obtains (3.32c). 

Next consider the case n > 1. The integrals in (D.3) are exactly the same as when n = 
except for the fact that the one in (D.3^) is followed by Z'^". Using the same arguments as 
before, we obtain 

Mi'o"'\z,t) = Q„-i(0/.XfcOZ-i +n'^;^^'\z,t)Z-^" + ... , (D.5a) 

H^X'^izJ) = "iQi{t)n'l^''\z,t)Z-' +n'iS'\z,t). (D.5^) 

The induction with (D.5), one can derive (D.2) for n>\. Sinnilarly, one obtains (D.2) for 
This completes the proof of (D.2). 

The above results imply that ju^^\z, t) = \ + 0{Z~^) as z ^ (oo, 0). In particular, computing 
the 0(7.~^) terms explicitly one obtains the first of (3.32fl). Similarly, using the same 
arguments, one can show that iJ^\z,t) = I + OiZ~^) as z — > (oo,0) and verify the second 
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of (3.32a). In the IVP, the integrals in the RHS of (D.3a) and (D.3fc) are absent, and the 
summation starts from n' - -oo. Hence in this case one simply obtains (3.17). 

The determination of the asymptotic behavior of /i„ (z,0 requires a shghtly different 
approach, since following the above steps for f/^\z,t), yields a 0(1) term involving the 
summation of Q„ in the RHS. To circumvent this difficulty, note that (A.5a) impUes /j^^ = 
(Z + Q„(0)"VfH_\ • For/}„(z,f) = C„i2'-^\z,t) we have 

/in -Z"'//„+i = -Q„/i„Z, (D.6) 

with finizj) —> I as «—> oo thanks to (3.25a) and (3.6). Introducing the auxiliary function 
%(z,0 = Z~"/i„(z,?), it is easy to check that *Pn(z,0 satisfies the equation *P„+i - ^n = 
ZZ~("+i'(Q„)1'„+i , which can be integrated to obtain the modified lost solution as 

Mz,t) = \-Z 2 Z''-''(Q„,_i(0/i„'(z,0). (D.7) 

n'=n+\ 

Then, applying the same Neumann series approach as described above to (D.7), one finds the 
asymptotic expansion for yuj,^^ as (3.18). Since i/^\z,t) is the same in the IVP and in the IB VP; 
this asymptotic behavior applies to both problems. 

Note that the above results also imply that a(z) and diz) are even functions in D±in and 
the following symmetries of M*: 

K,n(-Z't)^Kn(Z't), M;i2(-2,0 = -M:_i2(z,0, (D.8a) 

K,n(-z,t) = Kn(z,t), M-i2(-2,0 = -M-i2(-z,0. (D.8^) 



Appendix E. Independence of the solution on T 

The solution of a DDE does not depend on future values of the BCs. Hence, for any Tq < T 
the solution of the IBVP resulting from the RHP obtained by replacing T with To must be 
equivalent for all < ? < Tq to the solution of the IBVP obtained from the original RHP. We 
show next that is indeed the case because the RHP obtained from Tq and T are related. 

Let M„(z,f) satisfy the RHP (3.36), and let Mf "(z,f) and MJ™'(z,f) denote the restrictions 
of M„(z,t) to the domains D+m and D+om, respectively. Moreover, let A{z, To) and B{z, To) be 
the spectral coefficients obtained by replacing T with To in (E.3), and let ^^\z, t),..., ^^\z, t) 
denote the jump matrices obtained by replacing A{z,T) and B(z,T) with A(z,To) and 
B(z,To)- Finally, let M„(z, f) satisfy the RHP with the jump matrices J^\...,J^^ replaced 
by J^^', . . . , It is straightforward to see the relations 

m:'" = (I - ji'>) (I - ji'Y , M„-'" = M„-'° , 

M+out ^ M+out ^ m;-' = M-"' (I - 2fY (I - 4'^) • 

Now recall that ^„(f) can be obtained from the eigenfunctions via (3.40) or (3.32a) with 
7=1. Note also that /J^'^\z,t) enters M~™ via (3.35b). Below, we show that the matrices 
(I - J^'^) (I - J^' V and (I - J^^ V (I - J^^') are analytic and bounded for z g £>+in and z g £>+out, 
respectively. Since M„(z,?) = Mniz,t) for z g D-m, it then follows that the solutions qnit) 
obtained from M„ and M„ coincide. 

To show that (I - j),'^)(l - J^'^) ' is analytic and bounded for z e £>+in, note first that 

= ( I -^"^""'^"r (1/^%?-) -r*(i/z%ro)) j ^^^^ 

and the (1,2) component of (E.l) can be written as 

y r.^ 2n-2Mz)t A*(l/z\To)B*(l/z\T)-A*(l/z*,T)BWz*,To) 

A„(z) = vz e .(E.2) 

diz, T)d{z, To) 
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Now note that (3.28) and (3.31) define the scattering data A(z,r) and B{z,T) as 
Hence 

diz,T)d(z,To) 

Also, //q '^\z,f) satisfies the second column of the ?-part of the Lax pair (A.5) at n = 0: 

fiiiz, t) + 2ioj{z)ni{z, t) =Ho,n (z, t)iii{z, t) + J¥o,i2(z, t)n2{z, t) , (E.4a) 

fi2{z,i) = Ho,2iiz,t)fiiit,k)+Ho,22iz,t)fi2iz,t). (E.4fo) 

Then, introducing 

^i(z,0=/i;(l/z*,r)jui(z,0-//i(z,7^)/^;(l/z*,0e^'"^'^^^-^ (E.5a) 

niz, t) =fx*2il/z\ T)n2iz, t)-vfii (z, T)n\{l/z*, t) e2''^«(^-') , (E.Sb) 
we can rewrite the (1,2) component of (I - J^^')(l - J^'')~^ as 

2n„2iaj{z){To-t) 

X„iz) = ,, „, fi (z, To) . (E.6) 

d(z, T)d(z, To) 

It is therefore enough to show that (z, t) is analytic and bounded for z e D+ . The symmetries 
of Hoiz,t) [namely, Ho,i2(z,0 = vH*2i(l/z*,0 and Ho,ii(z,0 = Ho,22(l/z*,0] imply that 
{(f\,(p2y satisfies the f-part of the Lax pair (A.5b) with « = 0. Since ip\{z,T) = and 
Viiz, T)- 1, we then have the following linear integral equations 

ipi{z,t) = - /e2''^«(''-'>(Ho,ii.^i + Ho,i2^2)(z/)dr' , (E.7a) 

t 

T 

ip2(z,t) = 1 - / (Hoam + »o,22V2)(z,t')dt' , (E.7^) 
f 

From here one can show that (pi and (p2 are analytic and bounded for z e D+. As a 
result, the RHS of (E.2) is analytic and bounded for z e £>+. Thus (I - J^'^) (I - J^'V is 
analytic and bounded for z € D+in. The result for (I - J„ ) (I - J„ ) follows from symmetry 
considerations. 



Appendix F. Linearizable BCs for T koo 

Here we verify that (4.10) can be used to express r*(l/z*) also when T < oo. To do so, we 
use the same approach that we used to show that the solution of the IDNLS equation does not 
depend on T. Denote by Z„(z) the difl'erence between the contributions to the RHP obtained 
from T = oo and T < oo, namely: 

X„(z) = yz2"e-2'"(*(r(l/z*)-r:(l/z*)), (F.l) 

where r*(l/z*) is obtained by neglecting the second term in the RHS of (4.11). We can 
write (F.l) as 

X (z) = vz2«e-2''^(^)' Riz,T)-Roiz,T) 

d(z)do(z)/A*(l/z*,T)AUl/z*,T)' 

with Riz,T) = B*{llz\T)IA*{llz*,T) as before, and where Roiz) = fi;(l/z*)/A*(l/z*) is 
computed using only the first term in the RHS of (4.11) and doiz) = a(z)Al{l/z*,T) - 
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vb{z)Bl{llz*,T). Also, A;(1/z*,T) and Blil/z*,T) are defined by (4.9). Now, using (4.11), 
we find 

X„iz) = z2«e2'-(^X^-0/(l/z) l^iMl . (F.2) 

t/(z)A(l/z) 

In the solitonless case, however, we can assume that d(z) and A(l/z) never vanish in 5+in. 
Then the RHS of (F.2) is analytic and bounded in D+i„ due to the exponential term and now 
we know that the additional term in (4.1 1) does not affect the solution of the RHP. Note that 
/(1/z) has a pole at z = +1/;^^'''^. When 1 <x, or < -1, these points belong to D+m. Note, 
however, that since a(z) and b(z) are bounded in D+i„, if f(l/z) has a pole, A(l /z) does too, 
and hence the terms causing the poles in (F.2) to cancel out. 
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